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Abstract 

Lie group theory was originally created more than 100 years ago as a 
tool for solving ordinary and partial differential equations. In this article 
we review the results of a much more recent program: the use of Lie 
groups to study difference equations. We show that the mismatch between 
continuous symmetries and discrete equations can be resolved in at least 
two manners. One is to use generalized symmetries acting on solutions 
of difference equations, but leaving the lattice invariant. The other is 
to restrict to point symmetries, but to allow them to also transform the 
lattice. 
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1 Introduction 

The symmetry theory of differential equations is well understood, it goes back 
to the classical work of Sophus Lie and is reviewed in numerous modern books 
and articles [14, 21, 23, 34, 74, 104, 106-109, 194, 198, 208, 224-226, 236]. As a 
matter of fact, Lie group theory is the most general and useful tool we have 
for obtaining exact analytic solutions of large classes of differential equations, 
specially nonlinear ones. 

The application of Lie group theory to discrete equations is much more 
recent and a vigorous development of the theory only started in the 1990-ties 
[12, 17, 31, 42, 50-58, 63, 64, 79, 80, 86, 89-95, 112, 113, 118-120, 123-125, 127-130, 
134-145, 148-161, 166-171, 176, 189, 199-201 , 207, 209, 213-216, 218, 233, 237-239] . 

The purpose of this article is to provide a review of the progress made. 

In this whole field of research one uses group theory to do for difference 
equations what has been done for differential ones. This includes generating new 
solutions from old ones, identifying equations that can be transformed into each 
other, performing symmetry reduction, and identifying integrable equations. 
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When adapting the group theoretical approach from differential equations 
to difference ones, we must answer three basic questions: 

1. What do we mean by symmetries? 

2. How do we find the symmetries of a difference system? 

3. What do we do with the symmetries, once we know them? 

Let us first briefly review the situation for differential equations. 
Let us consider a general system of differential equations 

E a (x,u,u x ,u 2x ,...u nx ) =0, x G W,u G M. q ,a = 1, . . .N, (1.1) 

where u nx denotes all (partial) derivatives of u of order n. The numbers p, q, n 
and N are all nonnegative integers. 

We are interested in the symmetry group of system Ql.lfl . i.e. in the 
local Lie group of local point transformations taking solutions of eq. (jl.l|l into 
solutions of the same equation. Point transformations in the space X x U of 
independent and dependent variables have the form 

x = Aa(i, u), u = Vt\(x, u), (1.2) 

where A denotes the group parameters. Thus 

Aq(x,u)=x, tto(x,u)=u, 

and the inverse transformation (x,u) <— > (x,u) exists, at least locally. 

The transformations ljl.2|l of local coordinates in X x U also determine the 
transformations of functions u = f(x) and of derivatives of functions. A group 
of local point transformations of X x U will be a symmetry group of system 
<|1.1[1 if the fact that u{x) is a solution implies that u(x) is also a solution. 

How does one find the symmetry group (3? Instead of looking for "global" 
transformations as in eq. (|1.2(l one looks for infinitesimal ones, i.e. one looks for 
the Lie algebra g that corresponds to 0. A one-parameter group of infinitesimal 
point transformations will have the form 

Xi = Xi + \£i(x, u), u a — u a + \<j> a (x,u), (1.3) 
|A| < 1 1 < i <p, l<a<q. 

The search for the symmetry algebra 3 of a system of differential equations is 
best formulated in terms of vector fields acting on the space X x U of independent 
and dependent variables. Indeed, consider the vector field 

v 1 
X = y^^j{x,u)d Xi + ^2 4>a{x,u)d Ua , (1.4) 

i— 1 a—1 

where the coefficients £j and </> a are the same as in eq. (|1.3ll . If these func- 
tions are known, the vector field l|1.4fl can be integrated to obtain the finite 
transformations l|1.2|) . Indeed, all we have to do is to integrate the equations 

-jr = -JT = <l> a (x,u), (1.5) 
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subject to the initial conditions 

Xi |a=o= Xi u a \ x=0 =u a . (1.6) 

This provides us with a one-parameter group of local Lie point transformations 
of the form 11.2fl where A is the group parameter. 

The vector field i|1.4|) tells us how the variables x and u transform. We also 
need to know how derivatives like u x , u XXl . . . transform. This is given by the 
prolongation of the vector field X. 

We have 

I"' A ' V E{E<^< +E^ fc d UcWfc (1.7) 

ck i i,k 

+ £ ^ XkX 'du mi +■■■}, 
i,k,l J 

where the coefficients in the prolongation can be calculated recursively, using 
the total derivative operator, 

~\~ y>a,Xi&u a r ^ j a,x a Xi^u a :Ca ~t~ ^ j a,x a xtXi^u a ^ aJ:h ~\~ • • • (^-'^) 

(a summation over repeated indices is to be understood) . The recursive formulas 
are 

= D Xi (p a - (D Xi £ a )u a , Xa , <j) x a ' Xk = Dx k 0a - (D Xk £ a )u a , XiXa ,(1.9) 
^ XkXl = D XI <I%"> - (D Xl Z a )u a<XiXkXa , 

etc. 

The invariance condition for system (|l.lfl is expressed in terms of the oper- 
ator (|1.7|l as 

pr (B) ^ak=..^=o=0, a = l,...N, (1.10) 

where pr( n ) X is the prolongation (|1.7|) calculated up to order n (where n is the 
order of system (|1.1|) 1. 

Eq. (|1.10() provides a system of linear partial differential equations for the 
functions u) and <p a (x,u), in which the variables x and u figure as inde- 
pendent variables. By definition of point transformations the coefficients and 
4> a depend only on (x\, . . . , x p , Ui, . . . , u q ), not on any derivative of u a . The 
action of pr^X in eq. (|1.1U|) will, on the other hand, introduce terms in (|1.1U|) . 
involving the derivatives — fc ° - kp , k = k± + . . . + k p , 1 < k < n. We use equa- 

tions l|l. lfl to eliminate N (the number of equations) such derivatives. We then 
collect all linearly independent remaining expressions in the derivatives and set 
the coefficients of these expressions equal to zero. This provides the "determin- 
ing equations" : a set of linear partial differential equations for the functions 
£2(2, u) and (j) a (x,u). The order of the system of determining equations is the 
same as the order of the studied system ()1.1|) . however the determining system 
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is linear, even if the system Ijl.lfl is nonlinear. It is usually overdetermined and 
not difficult to solve. Computer programs using various symbolic languages ex- 
ist that derive the determining system and solve it, or at least partially solve 
it [35,88,221,224]. 

The solution of the determining system may be trivial, i.e. £j = 0, <j) a = 0. 
Then the symmetry approach is of no avail. Alternatively, the general solution 
may depend on a finite number K of integration constants. The Lie algebra of 
the symmetry group, the "symmetry algebra", for short, is then if -dimensional 
and must be identified as an abstract Lie algebra [87, 111, 181,205]. Finally the 
general solution of the determining equations may involve arbitrary functions 
and the symmetry algebra is infinite-dimensional. 

So far we have considered only point transformations, as in eq. I|1.2fl . hi 
which the new variables x and u depend only on the old ones, x and u. More 
general transformations are "contact transformations" , where x and u also de- 
pend on first derivatives of u [15,23,104,106,194,226]. A still more general class 
of transformations are generalized transformations, also called "Lie-Backhmd" 
transformations [15,114,194]. For these, the transformations are 

x = Ax(x,u,u x ,u xx ,...), u = fl\(x,u,u x ,u xx ,...), (1-11) 

involving derivatives up to an arbitrary, but finite order. The coefficients £j and 
(p a of the vector fields l|1.4fl will then also depend on derivatives of u. 

When studying generalized symmetries, and sometimes also point symme- 
tries, it is convenient to use a different formalism, namely that of evolutionary 
vector fields. Let us first consider the case of Lie point symmetries, i.e. vector 
fields of the form (|1.4|l and their prolongations (|1.7fl . To each vector field (|1.4fl 
we can associate its evolutionary counterpart X e , defined as 

X e = Q a {x,U,U x )d Ua , Q a = 4>a ~ £,jU a ,x ^ ■ (I- 12 ) 

The prolongation of the evolutionary vector field (|1 . 1 2|> is defined as 

w X e = Q,.i)„ ■ (>;, . • Q: t ' 0„„ • ... (1.13) 

The functions Q a are called the characteristics of the vector field. Observe that 
X e and pr X e do not act on the independent variables Xj . For Lie point symme- 
tries evolutionary and ordinary vector fields are entirely equivalent and it is easy 
to pass from one to the other. Indeed, eq. I|1.12f> gives the connection between 
the two. The symmetry algorithms for calculating the symmetry algebra g in 
terms of ordinary, or evolutionary vector fields, are also equivalent. Eq. I|1.10|l 
is simply replaced by 

pr W ^.k-.=B J ,=o,= 1 a=l,...N. (1.14) 

The reason that eq. (|1.10|) and i|1.14f> are equivalent is the following: 

pr< n > X e = w (n) X - ^Di. (1.15) 
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The total derivative Di is itself a generalized symmetry of eq. (|l.lf) . i.e., 



DiE a \ El=E2= ... =EN=0 ,= i = l,...p, a=l,...N. (1.16) 

Eqs. H1.15|l and l|1.16[) prove that systems (|1.1U|) and (|1.14|) are equivalent. 
Eq. (|1.16|) itself follows from the fact that DiE a = is a differential consequence 
of eq. i|l.l[l . hence every solution of eq. (|l.l|l is also a solution of eq. (|l.l(j|) . 

To find generalized symmetries of order k, we use eq. (I1.12[) but allow the 
characteristics Q a to depend on all derivatives of u up to order k. The prolonga- 
tion is calculated using eq. I|1.13|) . The symmetry algorithm is again eq. (|1.14|) . 

A very useful property of evolutionary symmetries is that Q a provide com- 
patible flows. This means that the system of equations 

< L17 > 

is compatible with system (| X . 1 |l . In particular, group invariant solutions, i.e., 
solutions invariant under a subgroup of are obtained as fixed points 

Q Q = 0. (1.18) 

If Q a is the characteristic of a point transformation, then ljl.181) is a system 
of quasilinear first order partial differential equations. They can be solved and 
their solutions can be substituted into (| 1 - 1 ft . yielding the invariant solutions 
explicitly. 

Many different extensions of Lie's original method of group invariant so- 
lutions exist. Among them we mention, first of all, conditional symmetries 
[22,67,68,147]. For differential equations, they were introduced under several 
different names [22,41,147,195] in order to obtain dimensional reductions of 
partial differential equations, beyond those obtained by using ordinary Lie sym- 
metries. 

Another valuable extension is the concept of partial symmetries. They corre- 
spond to the existence of a subset of solutions which, without necessarily being 
invariant, are mapped into each other by the transformation [39,40]. Further 
extensions are given by asymptotic symmetries [75,77], when extra symmetries 
are obtained in the asymptotic regime, or approximate symmetries [16,69,109] 
where one considers the symmetries of approximate solutions of a system de- 
pending on a small parameter. 

Let us now return to the problem at hand, namely symmetries of difference 
systems. We wish to study the continuous symmetries and use Lie algebra 
techniques. However, the equations are now discrete, i.e. they involve functions 
u(x) that are themselves continuous, but evaluated , or sampled, at discrete 
points. Several different approaches to this problem have been developed and 
are discussed below. 

In Section we discuss point symmetries of difference equations defined on 
fixed, nontransforming lattices [31, 79,86, 118-120, 123, 148, 151, 156, 167-170, 176, 
200,201,214,216]. The symmetry transformations are assumed to have the form 
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I|1.2[) . They must take solutions into solutions and the lattice into itself. This 
approach is fruitful mainly for differential-difference equations (DAE's), where 
not only the dependent variables, but also some of the independent ones are 
continuous. 

Section [3] is devoted to generalized point symmetries on fixed lattices [63, 
64,93, 128, 130, 139, 140, 145]. The concept of symmetry is generalized in that 
transformations act simultaneously at several points of the lattice, possibly in- 
finitely many ones. In the continuous limit they reduce to point transformations. 
This approach is fruitful mainly for linear equations, or equations that can be 
linearized by a transformation of variables. 

In Section 0] we consider generalized symmetries on fixed lattices [89-91,95, 
127] This approach is fruitful for discrete nonlinear integrable equations, i.e., 
nonlinear difference equations possessing a Lax pair. The symmetries are gen- 
eralized ones, treated in the evolutionary formalism. In a continuous limit they 
reduce to point and generalized symmetries of integrable differential equations. 

Point symmetries transforming solutions and lattices [17,50-52,54-58,141, 
143,157] are considered in Section [5] The transformations have the form (|1.2fl 
and they act on solutions and on lattices simultaneously. The lattices themselves 
are given by difference equations and their form is dictated by the symmetries. 
Their main application is to discretize given differential equations while preserv- 
ing their symmetries. 

A brief conclusion with an overview of the results presented and pending 
issues is given in Section [H] 

2 Point Symmetries of Difference Equations De- 
fined on Fixed, Nontransforming Lattices 

The essentially continuous techniques for finding Lie symmetries for differential 
equations can be extended in a natural way to the discrete case by acting just 
on the continuous variables [148, 151, 167-170,200], leaving the lattice invariant. 
Transformations of the lattice are considered only at the level of the group which 
itself is finite or discrete. Depending on the discrete equation we are considering 
we can have translations on the lattice by multiples of the lattice spacing and 
rotations through fixed angles, i.e., for example, by ^ when the partial difference 
equation is invariant with respect to the interchange of the two lattice variables. 

In Section IO we consider, as an example, the Lie point symmetries of the 
discrete time Toda lattice [127]. In Section we review the steps necessary 
to obtain continuous symmetries for differential-difference equations (DAE's) 
and then in Section [2. 31 we carry out the explicit calculations in the case of the 
Toda equation in 1+1 dimensions. In doing so we will present and compare 
the results contained in our articles on the subject, which correspond to both 
the intrinsic method and the differential equation method proposed earlier, and 
those of the global method introduced by Quispcl et al. [200] . 

Apart from the standard Toda lattice we will consider a Toda lattice equation 
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with variable coefficients and will show how, by analyzing its symmetry group, 
one can find the Lie point transformation which maps it into the standard Toda 
lattice equation with constant coefficients. In Section l2~4l we will present some 
results on the classification of nonlinear DAE's with nearest neighbour interac- 
tions. Finally, symmetries of a two-dimensional Toda lattice are discussed in 
Section 1231 

2.1 Lie point symmetries of the discrete time Toda equa- 
tion 

The discrete time Toda equation [100] is one of the integrable completely discrete 
partial differential equations (AAE) [11,99,101-103,131,186,187,190,196,217, 
219,220,227] and is given by 

A Toda = e «-.»-«».m+l _ e "».m +1 -»».m+» _ (2.1) 
- a 2( e «n-l, m+2 -tt„,m+l _ e «n,m+l-«»+l,m) = Q 

On the left hand side of eq. I|2.1|l we can easily obtain the second difference of 
the function u n ,m with respect to the discrete-time m. Thus, defining 

t = ma t \ v n (t)=u n , m ] ot = a\ (2.2) 

we find that eq. l|2.1ll reduces to the continuous-time Toda equation: 

4 2 L = V« - e""- 1 — + e^^ 1 = 0, (2.3) 

when er t — > and m — > oo in such a way that t remains finite. The Toda equa- 
tion (|2.3|l is probably the best known and most studied differential-difference 
equation. It plays, in the case of lattice equations, the same role as the Korteweg 
- de Vries equation for partial differential equations [230,231]. It was obtained 
by Toda [229] when trying to explain the Fermi, Pasta and Ulam results [62] on 
the numerical experiments on the equipartition of energy in a nonlinear lattice 
of interacting oscillators. As shown below in Section F4.2.2I ea. (|2.3|l reduces, 
in the continuum limit, to the potential Korteweg-de Vries equation. It can be 
encountered in many applications from solid state physics to DNA biology, from 
molecular chain dynamics to chemistry [206] . 

Eq. I|2.1|l can be obtained as the compatibility condition of the following over- 
determined pair of linear difference equations: 

1 gMn.m— «n,n+l 

^n-i,m + («+-- ae^.m+i-^,m (2.4) 

a a 

1pn,m+l = 1pn,m - ae Un ' m+1 ~ Un+1 - m tp n +l,m ( 2 -5) 

Let us consider the Lie point symmetries of the discrete-time Toda lattice l|2.1|) , 
on a fixed nontrasforming bidimensional lattice characterized by two lattice 
spacings in the two directions m and n, at and a x . 
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If the lattice is uniform and homogeneous in both variables, we can represent 
the lattice by the following set of equations: 

x n ,m = na x , t n>m = ma t . (2.6) 

Eqs. Q2.ti|) from now on will be denoted as ^-Lattice = 0. 

A Lie point symmetry is defined by giving its infinitesimal generators, i.e. 
the vector field 

Xn,m — £n,m {^n,m > ^n,m j ^n,m)^a: n . m ~t~ (2-7) 

and generates an infinitesimal transformation in the site (n, m) of its coordinates 
and of the function u(x niTnj t n ^ m ) = u n ^ rn . The action of 1|2.7|) on the discrete- 
time Toda lattice equation (|2.1|) is obtained by prolonging 1|2.7[) to all points of 
the lattice. The prolongation is obtained [143,148,151] by shifting (|2.7|) to these 
points 

prX = X n>m + X n +i jIn + X n>m +1 + Xn^n+2 + Xn-i^m+2- (2-8) 

The invariance condition then reads: 

prXA Toda \ {AToda=0ALattm=0) = 0, prXA Lat ti ce \(A Toda =o,ALatnoe=Q) = 0(2-9) 

The action of l|2.8|) on the lattice equation (|2.t)|) gives m = 0, r„ jm = 0. When 
we act with (|2.8fl on the Toda equation, we get 

e Un - m ~ Un - m+1 [^ n ,m - </>n,m+l] - e""' m + 1 " U '*'" l+2 [0„, m+ l - (j) n ,m+2\ ~ (2.10) 

-a a {e u »- I "»+ 3 - u -"»+ 1 [^ n _i, ro + a - (f> n , m +i] - 

_ e «n, m+1 -«» +I , m ^ n>m+1 _ m ]} = 

From the action of X on i|2.6[) the infinitesimal coefficients £ and t are zero 
and thus the variables x and t are invariant. If we differentiate eq. (|2.10|l twice 
with respect to u n ^ m +2 we get that 4> n , m = c\e Un - m + C2, where c\ and C2 are 
two integration constants (that can depend on x n ^ m and i„, m ). Introducing this 
result into eq. (|2.1U|I we get that ci must be equal to zero. Taking into account 
that, due to the form of the lattice, all points are independent we get that C2 
must be just a constant. 

To sum up, the discrete time Toda equation (|2.1fl considered on a fixed 
lattice has only a one-dimensional continuous symmetry group. It consists of the 
translation of the dependent variable u, i.e. u n ^ m = u n>m + A with A constant. 
This symmetry is obvious from the beginning as eq. i|2.1f> does not involve 
u n ,m itself but only differences between values of u at different points of the 
lattice. Other transformations that leave the lattice and solutions invariant will 
be discrete [138]. In this case they are simply translations of x and t by integer 
multiples of the lattice spacing a x and at ■ 
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The same conclusion holds in the general case of difference equations on 
fixed lattices. Lie algebra techniques will provide transformations of the contin- 
uous dependent variables only, through the transformations can depend on the 
discrete independent variables. 

In Section 0] we will see that the situation is completely different when gen- 
eralized symmetries are considered. In Section [S] we will consider transforming 
lattices, which greatly increases the role of point symmetries. 



2.2 Lie point symmetries of DAEs 

Let us now consider the more interesting case of DAEs. For notational simplic- 
ity, let us restrict ourselves to scalar DAEs for one real function u(n, i) depend- 
ing on one lattice variable n and one continuous real variable, t. Moreover, we 
will only be interested in DAEs containing up to second order derivatives, as 
those are the ones of particular interest in applications to dynamical systems. 
We write such equations as: 



Al 2 ) = Aft 



/. I 00 101 I 02 \ 

(t, n, u n+k | k=ag , u n+k ,t | k=ai , u n+k ,tt | k=a2 ) 



= 



aj < bj e z, (2.ii) 



with u n = u n (t). 

The lattice is uniform, time independent and fixed, the continuous variable 
t is the same at all points of the lattice. Thus to eq. 12.11f) we add 



(2.12) 



Examples of such equations which will be considered in the following are the 
Toda lattice equation (|2.3|) and the inhomogeneous Toda lattice [133] 



Al 2) = w M {n)-\w >t + \-^ 



-(n-lf + l 



w(n—l) — w(n) 



i 



e w(n)-w(n+X) = Q (2.13) 



We are interested in Lie point transformations which leave the solution set 
of eqs. I|2.11|l . (|2.12|) invariant. They have the form: 



f = A g (t,n,u n (t)), Un{t) = n g (t,n,u n (t)), 



n = n 



(2.14) 



where g represents a set of continuous or discrete group parameters. 

Continuous transformations of the form (|2.14|l are generated by a Lie algebra 
of vector fields of the form: 



X = T n (t,U n (t))d t + <pn(t,u n (t))d Un 



(2.15) 



where n is treated as a discrete variable and we have set n — n, when considering 
continuous transformations. 

Invariance of the condition (|2.12|) implies that r does not depend on n. 
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As in the case of purely differential equations, the following invariance con- 
dition 

P r( 2 )XA(?)| A c 2)=o = (2.16) 

must be true if X is to belong to the Lie symmetry algebra of A„ . The symbol 
pr^ 2 ) X denotes the second prolongation of the vector field X, i.e. 

n+b 

pr ( 2 ) X = T(t,u n )d t + ^2 (t>k(t,u k )d Uk 

k—n — a 
n+6i 

+ <#k(*>"*>«M) a <*M ( 2 - 17 ) 

k— n— ai 
n+b? 

k—n — a2 

with 

(j)l(t,u k ,u k ^ t ) = D t (f) k (t,u k ) - [D t T k (t,u k )]u ktt (2-18) 
(/) k (t,u k ,u ki t,u ki tt) = D t (j)\(t, u k ,u k ,t) ~ [D t Tk(t,Uk)]uk,tt (2-19) 

where D t is the total derivative with respect to t. 

Let us notice moreover that 0* and <j> u are the prolongation coefficients with 
respect to the continuous variable. The prolongation with respect to the discrete 
variable is reflected in the summation over fc. 

Eq. (|2.16(l is to be viewed as just one equation with n as a discrete variable; 
thus we have a finite algorithm for obtaining the determining equations, a usu- 
ally overdetermined system of linear partial differential equations for r(i,it n ) 
and (f)(n, t, w n j . 

A different approach consists of considering eq. (|2.11|) as a system of coupled 
differential equations for the functions u n {t). Thus, in general we have infinitely 
many equations for infinitely many functions. In this case the Ansatz for the 
vector field X would be: 

X = r(t, +E^(M"i(«)}i)S» k (i) (2-20) 

fc 

where by {uj(t)}j we mean the set of all Uj(t) and j and fc vary a priori over 
an infinite range. Calculating the second prolongation pr^ 2 ) I in a standard 
manner ( see eqs. (|1.7|l . 111. 90 1 and imposing 

pr^XA^\ A ^ =Q =0 Vn,j (2.21) 

we obtain, in general, an infinite system of determining equations for an infinite 
number of functions. Conceptually speaking, this second method, called the dif- 
ferential equation method in Ref. [151], may give rise to a larger symmetry group 
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than the intrinsic method we introduced before. In fact the intrinsic method 
yields purely point transformations, while the differential equation method can 
yield generalized symmetries with respect to the differences (but not the deriva- 
tives). In practice, it turns out that usually no higher order symmetries with 
respect to the discrete variable exist; then the two methods give the same result 
and the intrinsic method is simpler. 

A third approach [200, 201] consists of interpreting the variable n as a con- 
tinuous variable and consequently the DAE as a differential-delay equation. 
In such an approach u n+ k{t) = exp[^]{u„(i)} and consequently the DAE is 
interpreted as a partial differential equation of infinite order. In such a case 
formula (|2.1ti|) is meaningless as we are not able to construct the infinite order 
prolongation of a vector field X. The Lie symmetries are obtained by requiring 
that the solution set of the equation = 12. 3[) be invariant under the 
infinitesimal transformation 

t = t + €T n (t, u n (t)), 

h = n + eu n (t,u n (t)), (2.22) 

U n {t) = U n {t) + e(f) n (t, U n (t)) . 

To obtain conditional symmetries we add to eq. I|2.11|l a constraining equa- 
tion which we choose in such a way that it is automatically annihilated on its 
solution set by the prolongation of the vector field. Such an equation is the 
invariant surface condition 

A« = </> n (t,u n (t)) - T n (t,u n (t))u„ t t(t) = 0. (2.23) 

In general eqs. I|2.11|l and l|2.23[l may not be compatible; if they are, then 
eq. (|2.23() provides a reduction of the number of the independent variables by 
one. This is the essence of the symmetry reduction by conditional symmetries. 
Due to the fact that eq. (|2.23l) is written in terms of r and cf>, which are the 
coefficients of the vector field X, the determining equations are nonlinear. The 
vector fields do not form an algebra, nor even a vector space, since each vector 
field is adapted to a different condition [22,41, 147]. 

2.3 Symmetries of the Toda lattice 

Let us now apply the techniques introduced in Section l2~2*l to the case of eq. (|2.3f) . 
In this case eq. <|2.16|) reduces to an ovcrdctcrmincd system of determining 
equations obtained by equating to zero the coefficients of [tyt] fc , k = 0, 1, 2, 3 
and of v n ±±. They imply: 

t = at + d, <fi = b + 2an + ct, a, b, c, d real constants, (2.24) 

corresponding to a four dimensional Lie algebra generated by the vector fields: 

D = td t + 2nd Vn , f = d u W = td Vnl U = d Vn . (2.25) 
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The group transformation which will leave eq. (|2.3[l invariant is hence: 

v n (t) = v n (te^ 2 - A 3 ) + \ 2 {te~ x ^ 2 - A 3 ) + A 4 n + X 1 (2.26) 

where Xj, j = 1, 2, 3, 4, are real group parameters. 

To the transformation (|2.26l) we can add some discrete ones [136]: 

h = n + N NeZ (2.27) 

and 

(t, v n ) -» (-t, vn) ; (t, ««) -» (i, -«_„) • (2.28) 
We write the symmetry group of eq. Ij2.3|l as: 

# = 6 D i>(!l c (2.29) 

where €>n are the discrete transformations l|2.27l) . (|2.28|) and the invariant sub- 
group <&c corresponds to the transformation l|2.26ll . A complete classification 
of the one dimensional subgroups of <& can be easily obtained [151]. 
In fact, if we complement the Lie algebra l|2.25|l by the vector field 

2-£ (2.30) 

and require, at the end of the calculations, that the corresponding group pa- 
rameter be integer, the commutation relations become: 



[Z,D] = 2U; [D,T] = -T; [D,W] = W; [T,W]=U. (2.31) 

The one dimensional subalgebras are: 

{Z + aD + bU}, {Z + af + kW}, {Z + eW}, {Z + aU}, 

{f + cW}, {D + S}, {U}, {Z}, {W}, (2-32) 
{a, b, c) £ K; a ^ 0; k = 0, 1, -1; e = ±1. 

Nontrivial solutions, corresponding to reductions with respect to continuous 
subgroups & C $&Ci are obtained by considering invariance of the Toda lattice 
under {f + cW}, or {D + cU}. They are: 

1 " 

v n (t)=p--ct 2 -J2 l °s(Q-cj), (2.33) 

v n (t) =p + 2(n + c) log(t) - lo gb + (2c - l)j + 3% (2.34) 

j=o 

where p and q are arbitrary integration parameters. 

Reduction by the purely discrete subgroup, C 4I,d, given in eq. I|2.27|l 
implies the invariance of eq. I|2.3|l under discrete translation of n and makes it 
possible to impose the periodicity condition 

u(n + N,t) = u(n,t). (2.35) 
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This reduces the DAE (|2.3[l to an ordinary differential equation (or a finite 
system of equations). For example for N = 2, we get 

u tt = -4sinhu, (2.36) 

while for N = 3, we get the Tzitzeica differential equation [234,235]: 



Vtt 



(2.37) 



Let us now consider a subgroup (do C <& that is not contained in €>c, nor in 
#D! i- e - a nonsplitting subgroup of <©. A reduction corresponding to Z + aD + bU 
yields the equation 



F"(y) = e- b [exp(F(ye a ) - F(y) + a) - exp(F(y) - F(ye~ a ) 
where the symmetry variables F(y) and y are defined by: 
y = ter an , v n {t) = an 2 + bn + F(y). 
Using the subalgebra Z + aT + kW we get 



0] 



F"(y) 



o F(y+a)-F(y) 



o F(y)-F(y-a) _ * 



where 



y = t 



an, 



v n {t) = 



F(y)- 



(2.38) 
(2.39) 

(2.40) 
(2.41) 



Eq. H2.38fl can be called a "differential dilation" type equation; it involves 
one independent variable y, but the function F and its derivatives are evaluated 
at the point y and at the dilated points ye a and ye~ a . Eq. (|2.40(l is a differential 
delay equation which has interesting solutions, such as the soliton and periodic 
solutions of the Toda lattice (for k = 0). 

The other two nonsplitting subgroups give rise to linear delay equations 
which can be solved explicitely. 

This same calculation can also be carried out for the inhomogeneous Toda 
lattice (|2.13|) . The symmetry algebra is 



D = 2d f 



T = e"*/ 2 



dj- 



U 



(2.42) 



W = 2e' /2 a, 



u = a, 



These vector fields have the same commutation relations as those of the Toda 
lattice (12 .3|) . This is a necessary condition for the existence of a point transfor- 
mation between the two equations. In fact by comparing the two sets of vector 
fields, we get the following transformation which transforms a solution v n (t) of 
equation (j2.3|l into a solution w n (t) of eq. H2.13|) 



t = 2 log 



W n (t) = V n (t) - 21 71 - 



1 



2 

log(t) 



E 

3=0 



-(J-1) 2 + 1 



(2.43) 
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Let us now calculate the conditional symmetries of the Toda lattice. We 
assume that r is not zero and the determining equation reads: 



4>n,tt + 4>n,t4>n,v n + 2<t>n<t>n,tv n + 4>n[4> 2 n ,v n + 4>n4>n, v n v n )] 

+ [20„ - <t>n-l - 4>n+l}e" n ~ Vn + 1 + [4>n - <l>n-l][4>n,t ~ (/>n<f>n,v„ 



This implies 
with 



<f> n (t,v n (t)) =a(t)+p(t)n 



A« + P0,t = 0; a ,tt + /3a,t = °- 
Solving eqs. i|2.46[l we get 



<f>n(t, v n (t)) 



K + (^2^-+nK 3 j lanh 
K t + Ki + 2n 

t^T ' 



K ? 



(t - t ) 



, for K 3 ?0 



for Ka = 0. 



: 0. 

(2.44) 

(2.45) 
(2.46) 

(2.47) 



For K 3 = we get the results given in eq. (|2.24|) . For K 3 ^ an additional 
"symmetry" is given by <j> = K 3 n tanh(_K3t/2). This gives a new explicit solution 
of the Toda lattice: 



, (i) = uq + 2nK 3 log 



COhll( —^-t 



(2- 



2.4 Classification of differential equations on a lattice 

Group theoretical methods can also be used to classify equations according to 
their symmetry groups. This has been done in the case of partial differential 
equations [78] showing, for instance, that in the class of variable coefficient 
Korteweg-de Vries equations, the Korteweg-de Vries itself has the largest sym- 
metry group. The same kind of results can also be obtained in the case of 
differential-difference equations. 

Let us consider a class of equations involving nearest neighbour interactions 
[156] 

A„ =u nttt {t)-F n {t,u n -i{t),u n {t),u n+ i{t)) = 0, (2.49) 

where F n is nonlinear in Uk(t) and coupled, i.e. such that F n>Uh ^ for some 
k ^ n. 

We consider point symmetries only; the continuous transformations of the 
form (|2.14l) are again generated by a Lie algebra of vector fields of the form (|2.15|) . 
Taking into account the form of eq l|2.49|l . we have 



X = r(t)d t 



~T t + a n )u n + />„(/ ) 



du n (2.50) 
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with a n j — 0. The determining equations reduce to 



-TtttUn + b n>t t + (a« 



E 

fc=0,±l 



-T t )F n - rF n 



-r f + a n+k )u n+k + b n+k (t) 



0. (2.51) 



Our aim is to solve eq. 1)2. 51[) with respect to both the form of the nonlinear 
equation, i.e. F n , and the symmetry vector field X, i.e. (r(t), a n , b n (i)\ . In other 
words, for every nonlinear interaction F n we wish to find the corresponding 
maximal symmetry group d>. Associated with any symmetry group 41 there will 
be a whole class of nonlinear differential-difference equations related to each 
other by point transformations. To simplify the results, we will just look for the 
simplest element of a given class of nonlinear differential-difference equations, 
associated to a certain symmetry group. To do so we introduce so called allowed 
transformations, i.e. a set of transformations of the form 



t = i(t), n = n u n (t) =CL n (u n (t),t) 



(2.52) 



that transform equation l|2.49|l into a different one of the same type. By a 
straightforward calculation we find that the only allowed transformations l|2.52|l 
are given by 



t = t(t), n = n, u n (t) 



u n {t) + B n (t) 



(2.53) 



with B n (t), A n , t(t) arbitrary functions of their arguments. 

Under an allowed transformation equation (|2.49(l is transformed into 



with 



t(t) = F n (n,i,w n+ i(t),u n (f),w n _i(i)) 



F n (n,t, {ll n+ l(t), U n (t), U n -l(t)}) - B n ,u- 



(2.54) 



3§t 
4 rf 



lt_ut 
l .t 



A n u n {t) \ (2.55) 



and a symmetry generator H2.50JI into 
>(*)*,« , 1 



X = [r{t)t t ]d- t 



2 t. 



;*(*) 



A n 



-T{t)B n , t (t) + B n (t) ( -T t (t) +a n )+ b n {t) 



du n - (2.56) 
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We see that, up to an allowed transformation, every one-dimensional symme- 
try algebra associated to eq. (|2.49l) . can be represented by one of the following 
vector fields: 

X 1 =d t + a} n u n d Un X 2 = a 2 n u n d Un X 3 = b n (t)d Un (2.57) 

where a? n with j = 1, 2 are two arbitrary functions of n and b n (t) is an arbitrary 
function of n and t. The vector fields Xj, j = 1, 2, 3 are the symmetry vectors 
of the Lie point symmetries of the following nonlinear differential difference 
equations: 

Xi : u n> u = e a "*/n(£ n+ i,£„,£„_i), with = Uje' a i t 

(u-) a « 

X 2 ■ U n ,tt = Unfn(t,Tln+l,rin-l), with Tjj = — j- 

(UnP 

X 3 ■ Un,u = ^j^U n + f n (t, Cn+1, Cn-l) with Q = Ujb n (t) - U n bj(t). 

(2.58) 

These equations are still quite general, as they are written in terms of arbitrary 
functions depending on three continuous variables. More specific equations are 
obtained for larger symmetry groups [156]. 

The Toda equation (|2.3(l is included in a class of equations whose infinitesi- 
mal symmetry generators satisfy a four dimensional solvable symmetry algebra 
with a nonabelian nilradical. The interactions in this class are given by 

K(i,u„_i(t),M„(t),M„ + i(t)) = e" 2 T»£-T»/„(0 (2.59) 

where £ = (j n (t ) - 7„+ 1 (t ) )u n - 1 + (7„+i (t) - 7„_i (t) )u n + (7„_i (t) - j„ (t) ) u n+ i 
and the function j n (t) is such that 7„+i(<) ^= Jn(t) an d dl Qt = 0. The 
associated symmetry generators are: 

Xi = d Un , X 2 = d u X 3 = td Un , Y = td t + ln {t)d Un . (2.60) 

The Toda equation l|2.3|l is obtained by choosing j n (t) = 2n and f n (£,) = — 1 + 
ea?. Among the equations of the class (|2.49() . the Toda equation does not have 
the largest group of point symmetries. 

A complete list of all equations of the type (|2.49(l with nontrivial symmetry 
group is given in the original article [156] with the additional assumption that 
the interaction and the vector fields depend continuously on n. Here we just 
give two examples of interactions with symmetry groups with dimension seven. 
The first one is solvable, nonnilpotent and its Lie algebra is given by 

Xi = d Un , A 2 = (- l) n 8 Un , A 3 = td Un , 

Xi = (-l) n td u _ n , X 5 = (-l) n u n d Un , X 6 = d t , X 7 =td t + 2u n d Un . 

(2.61) 

The invariant equation is 

u„,tt = — • (2.62) 

U n -1 - U n+1 
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This algebra was not included in Ref [156] because of its nonanalytical depen- 
dence on n (in X 2 , and X 5 ). 

The second symmetry algebra is nonsolvable; it contains the simple Lie al- 
gebra sl(2,R) as a subalgebra. A basis of this algebra is 

X\ — d Un , X2 = td Un , X3 — b n d Un 

1 (2.63) 
X 4 = b n td Un , X 5 = d t , X 6 = td t + -u n d Un , X 7 = t 2 d t + tu n d Un 

with b nt t — 0, b n+ i ^ b n . The corresponding invariant nonlinear differential 
equation is: 

n ' tt [(b n +i ~ b n )u n _x + (6„_x - b n+1 )u n + (b n - 6„-i)u„+i] 3 

where 7„ and b n are arbitrary n-dependent constants. 

In [240] the integrability conditions for equations belonging to the class (|2.49|l 
has been considered. It has been shown that any equation of this class which 
has local generalized symmetries can be reduced by point transformations of 
the form 

u n = a n (t, u n ), t — 9(t) (2.65) 
to either the Toda equation Ij2.3|l or to the potential Toda equation 

Vf = e u "+ 1 - 2tl " +u "- 1 . (2.66) 

2.5 Symmetries of the Two Dimensional Toda equation 

Let us now apply the techniques introduced in Section ^. 2l to the two-dimensional 
Toda system (TDTS) 

&tdts = u n>xt - e""-i-"" + e""""^ 1 - (2.67) 

where u n is also a function of t and x. The TDTS was proposed and studied by 
Mikhailov [178] and Fordy and Gibbons [66] and is an integrable DAE, having a 
Lax pair, infinitely many conservation laws, Backlund transformations, soliton 
solutions, and all the usual attributes of integrability [1,8,33]. 

The continuous symmetries for eq. (|2.67|) are obtained by considering the 
infinitesimal symmetry generator 

X = £ n (x,t,u n )d x + r n (x,t,u n )d t + 4> n (x,t,u n )d Un . (2.68) 

From the determining equation (|2.16(l we get 

T n = f(t), = h{x). <t> n = [h, x + f, t ) n + g(t) + k{x), (2.69) 

where f(t), g[t), h{x) and k(x) are arbitrary C°° functions. A basis for the 
symmetry algebra is given by 

T(f) = ffldt + nf tt d Un , U(g) = g(t)d Un , (2.70) 
X(h) = h(x)d x + nh x d Un , W{k) = k{x)d Un , 
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where, to avoid redundancy, we must impose k tX ^ 0. 
The nonzero commutation relations are 

\T{h\T{h)\=T{hh,t~h,th\ [T(f),U(9)] = U(fg, t ), (2.71) 
[X(ht), X(h 2 )} = X{hth 2 , x - h ltX h 2 ), 



cU(l), (hk, x ) tX = Q, hk, x = c. 

Thus {T(f), U (g)} form a Kac-Moody-Virasoro u(l) algebra, as do {X(h), 
W(k),h^ x 0, U(l)}. However the two u(l) algebras are not disjoint. This 
Kac-Moody-Virasoro character of the symmetry algebra is found also in the 
case of a (2 + 1) -dimensional Volterra equation. It is also characteristic of 
many other integrable equations involving three continuous variables, such as 
the Kadomtsev-Petviashvili or three-wave equations [36,45,46,146,177,197]. 
From the symmetry algebra we can construct the group of symmetry trans- 
formations which leave the TDTS l|2.67|l invariant and transform noninvariant 
solutions into new solutions. Moreover, we can use the subgroups to reduce the 
TDTS (|2.67|) to equations in a lower dimensional space. 

By adding to l|2.67|) the equation 

Mi ) = €n(%,t,u n )u niX +T n (x,t ) u n )u n ,t ~ </> n (x,t,Un) = (2.72) 

we obtain the conditional symmetries. It is easy to show that conditional sym- 
metries in the intrinsic method do not provide any further symmetry reduction. 
It is, however, worthwhile to notice that conditional symmetries in the differen- 
tial equation method, when 

X = £ n (x,t,Uj(x,t))d x + T n {x,t,Uj{x,t))d t + J~] <f>k(x, t, Uj(x, t))d Uh , (2.73) 

fc 

contain the Backlund transformation of the TDTS. In fact, by choosing £„ = 0, 
T n = 1, eq. (|2.73(l reduces to An = u nt t — <p n {x, t, Uj) — and the determining 
equations are solved by putting 

(f> n (x,t, Uj ) = / rM (M) +a[e tl "" / " +1 - e 11 "- 1 "/"] (2.74) 

where a is a real constant and /„(ir,£) = u n is a set of functions which solves 
the TDTS (|2~F7|) . The Backlund transformation for the TDTS can be indeed 
written as [66] 

u n , t - Un,t = a[e""+! - e"- 1 " 8 "], (2.75) 

^"n— l.x &n.x [^ ^ ]■ 

a 

3 Generalized Point Symmetries on a Fixed Uni- 
form Lattice 

We saw in Section [5] that point symmetries for purely difference equations on a 
fixed uniform lattice do not provide very powerful tools, though they work well 
for differential - difference equations. 
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Here we shall consider an approach that makes use of a certain type of 
generalized symmetries that act simultaneously at more than one point of the 
lattice. We call them "generalized point symmetries" , because in the continuous 
limit they reduce to point symmetries. 

This approach is directly applicable to linear difference equations, indirectly 
to nonlinear equations that are linearizable by a change of variables. 

The underlying formalism is called " umbral calculus" , or " finite operator 
calculus". Its modern development is mainly due to G.C. Rota and his col- 
laborators [210-212]. For a review article with an extensive list of references, 
see [47]. Umbral calculus has been implicitly used in mathematical physics 
[63,64,82,128,130,145]. The only explicit use in physics that we know of is 
in [37,48,139,140]. 

3.1 Basic Concepts of Umbral Calculus 

Definition 1. A shift operator T„ is a linear operator acting on polynomials 
or formal power series f{x) in the following manner 

Taf(x) = f(x + a), a e K. (3.1) 

For functions of several variables we introduce shift operators in the same 
manner 

T(7if • ■ - 1: 3^+1 • ■ - %n) 

= fi^l: • ■ • ! "i" > ■ • ■ ; (3-2) 

In this section we restrict the exposition to the case of one real variable 
x E K. The extension to n variables and other fields is obvious. We will 
sometimes drop the subscript on the shift operator T when that does not give 
rise to misinterpretations. 

Definition 2. An operator U is called a delta operator if it satisfies the 
following properties, 



1) It is shift invariant; 

2) 

3) 



T a U = UT a , Vct e R, 
Ux = c^0, c = const, 
Ua = 0, a = const, , 



(3.3) 
(3.4) 
(3.5) 



and the kernel of U consists precisely of constants. 
Important properties of delta operators are: 

1. For every delta operator U there exists a unique series of basic polynomials 
{P n (x)} satisfying 

P (x) = l P„(0)=0, n>l, UP n (x)=nP n -i(x). (3.6) 
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2. For every delta operator U there exists a conjugate operator (3, such that 

[U,x/3] = 1. (3.7) 

The operator (3 satisfies 

f3={U')- 1 , U' = [U,x]. (3.8) 

Eq. I|3.7|l can be interpreted as the Heisenberg relation between the delta oper- 
ator U and its conjugate x(3. 

We shall make use of two types of delta operators. The first is the ordinary 
(continuous) derivative operator, for which we have 

U = d x , 13=1, P n {x)=x n . (3.9) 

The second is a general difference operator which we define as 

1 m 

U = A = -V a k T*, l,meZ, Km, (3.10) 
a z — ' 

k=l 

where a k and a are real constants and T a is a shift operator. In order for A to 
be a delta operator, we must impose 

m 

J2^k = 0. (3.11) 

We shall also require that the continuous limit of A be d x . This imposes a 
further condition, namely 

TCI 

S £ j ka h = \. (3.12) 

k=l 

Than eqs. (ESJ, . • • , (E3 hold, with c = 1. 

If eqs. H3.11JI . I|3.12(l are satisfied, we shall call A a difference operator of 
order m — I. Acting with A on an arbitrary smooth function f{x) we have 

A f( x ) = -Y j a k f{x + ka) = -V^— pV 7 n, (3.13) 

fc=i «=0 

m 

In = ^a fc fc", 7 = 0,7i = l. 



k=l 



Thus we have 



ni 

n=2 



For m - ! > 2 we can impose further conditions, namely 

7„ = 0, n = 2, 3, . . . , m - I, (3.15) 
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and obtain 

A= jL + 0(a m - 1 ). (3.16) 
ax 

From i|3.8|) the operator f3 conjugate to A is 

m 

(3 = (Y^aukT*)- 1 (3.17) 

k=l 

and the basic polynomials are 

P n (x) = [x/3) n l. (3.18) 

It was shown in [139] that the ea. (|3.18|l yields a well defined polynomial of order 
n in x. For any differential operator A all coefficients in P n {x) are finite and 
each involves a finite number of positive powers of shifts in a. 

The simplest examples of difference operators and the related quantities are 

1. 

A+ = — — , (3 = T-\ (3.19) 
a 

Pn{x) — (x) n = x(x — a)(x — 2a) . . . (x — (n — 1)<t), n > 1 
Order m — l= 1 . 

2. 

1 - T^ 1 

A" = , = T, (3.20) 

a 

P n {x) = (x) n = x(x + a)(x + 2a) . . . (x + (n - 1)<t), n>l 
Order m — I = 1 . 

3. 

Pin{x) = x 2 (x 2 - 4a 2 ) ... (x 2 - (2n - 2) 2 a 2 ), n > I, 
P2n+i(x) = x(x 2 -a 2 )... (x 2 - [2n - 1) V), n > 1. 

Order m — I = 2. 

For any A we have Pq — 1 and Pi = x. 

An important tool in the umbral calculus is the umbral correspondence. This 
is a bijective mapping M between two delta operators U\ and L^- This will 
induce a mapping between the corresponding operators f3\ and /?2, and also 
between the corresponding basic polynomials: 

U, A U 2 , fa A /3 2 , P^) A if). (3.22) 
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Let us now consider linear operators L(xf3, U) that are polynomials, or for- 
mal power series in x/3 and U. Since the umbral correspondence preserves the 
Heisenberg commutation relation (|3.7|) . it will also preserve commutation rela- 
tions between the operators L(x(3, U). In particular, we can take 

U 1= d x , fo = l, U 2 = A, fo = P, (3.23) 

where A is any one of the difference operators <|3.1U[) and (3 is as in (|3.17|l . 
Let Ai be a Lie algebra of vector fields of the form 

n 

X a = y^a° a (x 1 , ■ ■ ■ ,x n )d Xj , [X a ,Xp] = c i a pX 1 . (3.24) 
i=i 

The umbral correspondence will map this algebra onto an isomorphic Lie algebra 
A2 of difference operators 

n 

K = J2<( X ^-'- ^n(3n)A Xj . (3.25) 
i=i 

3.2 Umbral Calculus and Symmetries of Linear Difference 
Equations 

Lie point symmetries of linear differential equations can be expressed in terms of 
commuting linear operators. Indeed, let us consider a linear differential equation 

Lu = (3.26) 

where L is some linear differential operator. The Lie point symmetry algebra 
of this equation can be realized by evolutionary vector fields of the form (|1.12|) , 
satisfying eq. I|1.14fl . If (|3.2ti|) is an ordinary differential equation of order 
3 or higher, or a partial differential equation of order 2 or higher, then the 
characteristic Q of the vector field i|1.12fl will have a specific form 

Q = f(~x) + Xu, (3.27) 

where /( x ) is the general solution of eq. (|3.26l) and X is a first order linear 
operator of the form [20] 

n 

X = '%2$t(2)d xt -4>&), (3.28) 
i=i 

satisfying 

[L,X]u\ Lu=0 = 0. (3.29) 

Using the umbral correspondence, we can carry this result over to a class 
of linear difference equations. Indeed, let us associate a difference operator Lu 
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with L by the umbral correspondence d Xi — > A I; , a;, — > x$i. Eq. I|3.2(j[l is 
replaced by a difference equation 

L D u = 0. (3.30) 

The analogue of eq. I|3.29|) will hold, namely 

[L D} X d ]u\ LdU=o = 0. (3.31) 

The difference operator Xjj will not generate point transformations taking so- 
lutions into solutions (they cannot be integrated to point transformation in 
general). They do, however, provide commuting flows, i.e. difference equations 
compatible with eq. I|3.30|l . 

3.3 Symmetries of the Time Dependent Schrodinger Equa- 
tion on a Lattice 

As an example let us apply the umbral correspondence to the free time - depen- 
dent Schrodinger equation. We obtain the difference equation 

1 " 

L D ^ = 0, L D = iA t + ~J2&x k x k - (3-32) 

k=l 

The symmetries of eq. <|3.32[1 will be represented in terms of the "discrete" 
evolutionary vector fields 

XE = Q D d^ + Q* D d^, Q D = r)D-T D A t ip-&A Xk i>, (3.33) 

where rjn 1 td and ^ D are functions of tj3t, XjPj and ip, ip* (the * indicates 
complex conjugation). 

As in the continuous case, it can be shown [139], that the characteristic Q 
will in this case have the form 

Qd = x(x k (3k, tfy) + iX D i>, (3.34) 

where x is the general solution of eq. H3.32JI . The first order difference operator 
Xd satisfies 

n 

X D = i[r D {tf3 t )A t - i VD ], (3.35) 

fe=i 

[L D ,X D }^\ LD i ) =o = (3.36) 

Explicitating and solving eq. (|3.36(l . we obtain a difference realization of the 
Schrodinger algebra, first obtained in the continuous case by Niederer [183]. A 
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basis for this algebra is given by the following operators: 

P =A t , Pj = A Xj , L j , k = (x j /3 j )A Xk -(x k l3 k )A Xj , (3.37) 

% " 1 

B k = (t(3 t )A Xk - -x k (3 k , D = 2(t/3t)A t - ^{x k f3 k )A Xk + -, 

k=l 

n . n 

C = (tp t ) 2 A t + Y^Wt){x k f3 k )A Xk + -tfr ~ j 5I(x fc /3 fc ) 2 , W = i. 

k=l k=l 

Comparing with the continuous limit (or using the umbral correspondence) , 
we see that Pq, Pj correspond to time and space translation, Lj tk to rotations, 
B k to Galilei boosts, D and C to dilations and projective transformations and 
W to changes of phase of the wave function. 

3.4 Symmetries of the Discrete Heat Equation 

As a further example let us consider the discrete heat equation 

A xx u(x,t) - A t u(x,t) = 0. (3.38) 

Eq. H3.38f) is a linear partial difference equation on a two dimensional lattice. 
Floreanini, Negro, Nieto and Vinet showed in [63] that (|3.38(l has a group of 
symmetries isomorphous to that of the continuous heat equation 

u tXX (x,t) -u jt {x, t) = 0. (3.39) 

This result can be easily recovered using the umbral calculus [128, 145]. 

Since (|3.38|) is linear, the symmetries are obtained by considering an evolu- 
tionary vector field of the form 

X e = Qd u = (rA t + ^A, + f)ud u (3.40) 

and the determining equation is give by 

A t Q - A xx Q\ A ^ u = AtU = Q (3.41) 

whose explicit expression is 

A t {iA x u + rA t u + fu) - [A xx (£A x u + rA t u + fu)\ Axxll=AtU = 0. (3.42) 

Defining D x f = [A, /] it is immediate to see that (D x f)l = (A/)l and for the 
operator D the Leibnitz's rule is satisfied: D x fg — (D x f)g + f(D x g). In term of 
the operator D we can split eq. (|3.42|l into the following overdetermined system 
of equations 

D x t = 0, D t r - 2D X £ = , . 

AC - D xx i - 2D x f = 0, D t f - D xx f = 0. { ^ 6) 
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From the umbral correspondence the solution of eq. (|3.43|) is 



t = T 2 (tp t f + ntfc + r 

£ = \{ti + 2T 2 tl3 t )x(3 x + ^tfc + h 



(3.44) 



where tq, t\, t 2 , £i, £o an d fo a- re arbitrary constants, functions of the lattice 
spacing and shifts. By a suitable choice of these constants, we get the following 
representation of the symmetries 



For any choice of the A and the corresponding [3 operators we have a different 
representation of the symmetry algebra of the heat equation. 

4 Generalized Symmetries on Fixed Uniform Lat- 
tices 

4.1 Generalized symmetries of difference equations 

We consider here the construction of generalized symmetries for difference equa- 
tions and show the structure of the infinite dimensional Lie algebra of point and 
higher symmetries. 

Generalized symmetries are symmetries whose infinitesimal generators de- 
pend on more than one point of the lattice and on derivatives with respect to 
the variables which vary in a continuous way. E. Noether was the first to notice 
in 1918 [191] that one can extend symmetries by including higher derivatives of 
the dependent variables in the transformation. They are more rare than point 
symmetries. However, we can obtain an infinite number of generalized symme- 
tries [194] when the system is integrable [1,8,33,59,182,192,241], i.e when the 
equations can be written as the compatibility condition for an overdetermincd 
system of linear equations (the Lax pair) and can be linearized either directly 
or via an Inverse Scattering Transform. 

Lie symmetries, both point and generalized [194], can be obtained as flows 
commuting with the equation at study. We introduce here just the minimal 
notions on integrability necessary to be able to construct the commuting flows. 
For more on integrability see [1,8,33,59,182,192,230,231,241]. 

A nonlinear partial differential equation for one function of two variables 



P = (A t u)d u , Pi = (A x u)d u 

W = ud u , B = [2tj3 t A x u + x[3 x u]d u 

D = [2t[i t A t u + xf3 x A x u + ^u]d u 

K = [(tp t ) 2 A t u + t xf3 x A x u + {^(xp x ) 2 + hp t )u]d, 



(3.45) 



E(x,t,u(x,t),u x (x,t),u t (x,t), . . .) = 



(4.1) 
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is said to be integrable if it can be characterized by a Lax pair [121] 

L(u)ip = \if>, (4.2) 
^ t = -M{u)i>, (4.3) 

a system of linear equations compatible only on the solution set of solutions of 
eq. Q4.1|) . In eqs. Q4.21 I4.3[l A is an eigenvalue, L and M are two linear operators 
in x with coefficients depending on u but not on A. The function ip, often called 
the wave or spectral function, depends on the independent variables (x,t) and 
on A. If X t = than A is an integral of motion, together with all the functions 
which depend only on it. In the particular case when eq. 14.21) represents the 
Schrodinger Spectral Problem 

ip,xx + [A - u]ip = 0, A = fc 2 , (4.4) 

for u(x) vanishing at infinity, the solution of eq. I|4.4|) has the following asymp- 
totic behaviour: 

V>(x,fc)-> b(k)e lkx + a{k)e~ lkx , (x -► +oo), (4.5) 

i/}(x,k)—* e~ ikx , ^ x _^ _ 0O )_ (45) 

It can be easily proved that, when u(x, t) evolves according to the KdV equation 
[33], the function a(k) is conserved, i.e. a(fc) = 0. Thus, an infinite set of 
conservation laws can be obtained by expanding the transmission coefficient 
l/a(A) in powers of A. 

The compatibility of eas. H4.2l 14. 3|) implies the Lax equation 

L t (u) = [L(u),M(u)] (4.7) 

if A t = or 

L t (u) = [L(u), M(u)\ + f(L(u), t) (4.8) 

for A t = f(X,t). Here f(z,t) is an entire function of its first argument. 

Eq. I|4.7|) or eq. (|4.8|) is equivalent to the nonlinear partial differential 
equation (|4.ip . In the first case we say that the equation l|4.1l) is obtained as 
an isospectral deformation of l|4.2|l while in the second case by a non-isospectral 
deformation. For isospectral deformations the existence of a Lax pair implies 
the existence of an infinite number of conserved quantities while in the non- 
isospectral case these quantities will have a prescribed evolution in time. 

The solution of the Cauchy problem for these systems is obtained by solving 
the Inverse Scattering Problem [1,8,33]. The integrability of the isospectral 
systems has also been confirmed by the Painleve test, by deriving the global 
properties of its solutions through singularity analysis [117]. There is strong 
evidence, supported by many results, that all integrable equations have the 
Painleve property, i.e. all solutions are single- valued around all movable sin- 
gularities. In the case of ordinary differential equations this idea was used by 
Kowalevski [115,116] to identify integrable rotating tops. For partial differ- 
ential equations it lead to the Ablowitz, Ramani and Segur conjecture [6, 7] 
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which states that any ordinary differential equation which arises as a (symme- 
try) reduction of an integrable partial differential equation possesses the Painleve 
property, possibly after a transformation of variables. This conjecture has never 
been proven but there are no known counterexamples. 

Each integrable nonlinear partial differential equation 14. be it isospec- 
tral, or non-isospectral, is characterized by an operator M(u), defined up to 
an arbitrary constant matrix. To each linear problem (|4.2|l we can associate a 
recursion operator C and thus a hierarchy of nonlinear partial differential equa- 
tions [33]. The operators M{u) appear in a sequence constructed in terms of L. 
The construction of C from the linear problem 1)4. 2f> is algorithmic [3,28,29,32]. 
The solution of the Cauchy problem for the integrable equation l|4.1|l is obtained 
by performing the Spectral Transform and solving the Inverse Spectral Problem 
for the linear problem l|4.2l) . Starting from a given initial potential u(x,t — 0) 
we can characterize the spectrum S[u] of eq. 1)4. 2J) . Solving the linear equation 
for the spectrum we can obtain the potential u(x,t) for any t > by solving 
the inverse problem. Each equation of the hierarchy 

Ej{{x,t,u{x,t),u x {x,t),ut(x,t), . . .) = (4.9) 

is characterized by a different linear evolution of its spectrum S^it)]. At the 
level of the spectrum, it is easy to prove that from the hierarchy of equations as- 
sociated to a given linear equation l|4.2|l we can obtain a set of commuting flows. 
Thus every equation of the hierarchy (|4.9|l can provide a set of evolutionary, in 
general generalized, symmetries [65,70-73]. We have described here this proce- 
dure in the case of partial differential equations, where it was first introduced. 
This procedure has then been extended to the case of differential-difference and 
difference-difference equations [4, 5, 12, 38, 49, 60, 61, 166, 175, 185, 188, 222, 223]. 
In the case of a differential-difference equation 

E{n,t,u n {t),u n+1 (t),u n>t (t),...) = (4.10) 

the linear operators L and M that characterize it, depend on the shift operators 
in the discrete variable n. The Lax equations 1)4.71 14.8fl are still valid. The 
recursion operator C and the Lax pair will depend on the shift operators, instead 
of x derivatives. This implies that the higher equations of the hierarchy 

Ej(n,t,u n (t),u n+1 (t),u nit {t), . . .) = (4.11) 

and the higher symmetries will depend on points further away from the point 
n instead of depending on higher derivatives. The situation is slightly different 
in the case of difference-difference equations 

E(n, m, U„ tm , u„+i, m , U n>m+ i, ...) =0. (4-12) 

The Lax pair in this case involves two linear functions L„.,„ and M„ )m of the 
shift operator in n with coefficients depending on u n ^ m The the linear equation 
l|4.3[) is replaced by 

VVi,m+l = -M„, m (tt)V'n,m- ( 4 -13) 
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The Lax equation in the isospectral regime (X rn +i = A m ) now reads: 

L n ^ m +iM nm = M nm L nm (4-14) 

and in the nonisospectral case, when A m +i = /(A m ), with /(.) an entire function 
of its argument, 

L n ,m+l^n,m = Mn : mf{L n ,m)- (4-15) 

Few results are known on generalized symmetries of difference-difference equa- 
tions [127,172,184]. In [127], using the techniques presented below, we computed 
the point and generalized symmetries for the discrete time Toda Lattice (|2.1|) . 
The symmetries are provided by flows commuting with it. We are considering 
Lie symmetries, so these flows are given by differential difference equations. In- 
deed the group parameter is a continuous variable. So the symmetries of the 
discrete time Toda Lattice (12. If are given by the Toda Lattice hierarchy of non- 
linear differential-difference equations l|4.20|l , with the evolution not in the time 
variable, but in the group parameter. 

The Painleve property for integrable PDEs is not directly applicable. Several 
approaches have been proposed in the recent years [2,18,43,44,83-85,96-98, 
193,202-204,228,232]. The first is the singularity confinement proposed by 
Grammaticos et. al. [84]. iln this approach one requires that any singularity of 
the solution of a discrete system be confined and should not propagate. This 
criterion is not sufficient, a system can satisfy the singularity confinement, but 
give rise to numerical chaos [97]. More refined techniques have been proposed 
later. Among them let us mention that of Conte and Musette [43] in which 
they propose a natural extension of the Painleve property to discrete systems, 
that of Viallet and Hietarinta [18, 96-98] based on the study of growth and 
complexity properties and that of Ablowitz et. al. [2] based on the study of 
the asymptotic behaviour using the Nevanlinna theory of value distribution and 
growth of meromorphic functions. 

In the following we will construct an infinite class of symmetries for the Toda 
lattice. We will construct formally all the generalized symmetries and present 
explicitly the simplest examples. We will show the structure of the algebra 
of the generalized and point symmetries, using the one-to-one correspondence 
between the configuration space and the spectral space, where all equations 
are linear. In the continuous limit the Lie algebra of the point symmetries 
of the potential KdV l|4.70[) are recovered, as are the generalized symmetries. 
Finally, we present the Backlund transformations for the Toda lattice and show 
their relation to the symmetries. For the corresponding results for the Volterra 
equation, the completely discrete Toda lattice and for the Discrete Nonlinear 
Schrodinger equation we refer to the literature [38,89,90,95, 127, 132]. 

4.2 The Toda system and its symmetries 

The Toda equation (t2~31l can be rewritten in the form of a system 

a n (t) = a n (t)(b n (t)-b n+1 (t)), b n (t)=a n - 1 (t)-a n (t), (4.16) 
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where 

b n = v n , a n = e v -- v '>+ 1 . (4.17) 
It can be associated to the discrete Schrodinger spectral problem 

tp(n - 1, t; A) + b n ip(n, t; A) + a n ij)(n + 1, t; A) = Xtp(n, t; A). (4.18) 

The time evolution of the wave function ip(n, t; A) is given by 

i> t (n, t; A) = -a„ (t)ij(n + 1, t; A). (4.19) 

For the point symmetries of the Toda equation see Section |2~51 
To the spectral problem (|4.18|) we can associate a set of nonlinear differential 
difference equations (the Toda Lattice hierarchy) 

where f\{C, t) is an entire function of its first argument and the recursion oper- 
ator £, is given by 

C { Pn \ = { Pnbn+1 + a «(<7» + 1n+l) + (b n - b n+ l)s n 

The operator (|4.21|) was firstly obtained by Dodd [49] and by Bruschi, Levi and 
Ragnisco [26] . In eq. (|4.21|) s n is a solution of the nonhomogeneous first order 
equation 

S„+l = -^—(s n -Pn)- (4.22) 

For any equation of the hierarchy 14.20J1 we can write down an explicit evolution 
equation for the function ip(n, t; A) [27,28,30] such that A does not evolve in time. 
This is possible if the following boundary conditions 

lim a n — 1 = lim b n = lim s n = 0, (4.23) 

| TO | — >oo |n| — >oc |n| — >oo 

are imposed on the fields a n , b n and s n . The boundedness of the solutions of 
eas. (14.2311 was not required in the literature [25,29,30] but it is necessary to 
get a hierarchy of nonlinear differential-difference equations with well defined 
evolution of the spectra. We can than associate to ea. (|4.18|> a spectrum S[u] 
[27,60,61,175,230] defined in the complex plane of the variable z (X — z + z" 1 ): 

S[u] = {R(z,t), zeC^z^Cjit), \ Zj \<l,j = l,2,...,N}, (4.24) 

where R(z,t) is the reflection coefficient, Ci is the unit circle in the complex 
z plane, Zj are isolated points inside the unit disk and Cj are some complex 
functions of t related to the residues of R(z, t) at the poles Zj. When a n , b n and 
s n satisfy the boundary conditions 14.23fl . the spectral data define the poten- 
tials in a unique way. Thus, there is a one-to-one correspondence between the 
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evolution of the potentials (a n ,b n ) of the discrete Schrodinger spectral problem 
H4.18fl . given by eq. I|4.20|) and that of the reflection coefficient R(z, t), given by 

^M =/i/l (A,t)^,t), v = z- y -z. (4.25) 

In eq. H4.25fl an d below, 4z denotes the total derivative with respect to t. 

The Toda system is obtained from eq. (|4.20|) by choosing fi(X,t) = 1 and 
thus the evolution equation of the reflection coefficient is given by 

^M=M«(*,t). (4-26) 

The symmetries for the Toda system 1)4.16(1 . or the Toda equation i12.;}[l are 
provided by all flows commuting with the equation itself. Let us introduce the 
following denumerable set of equations 



£fc ~n^n "„+l, (4 27) 

bn,e k ) \ ln-1 — a n J 

Here k is any positive integer and e& is a variable. We can associate to the 
equation (|4.27|) an evolution of the reflection coefficient 

^ = fiX k R. (4.28) 

It can be verified directly that the flows l|4.26[) and (|4.28|) commute and hence 
the same must be true for the corresponding equations, i.e. for the Toda system 
l|4.16[) and the equations l|4.27ll . This implies that eqs. 14.27fl are symmetries of 
the Toda system and consequently tk is a group parameter. From the point of 
view of the spectral problem l|4.18|l . eqs. I|4.27[) correspond to isospectral defor- 
mations, i.e. we have A efc = 0. For any ek, the solution of the Cauchy problem 
for eq. (|4.27|) . provides a solution of the Toda system (|4.16|l [a„(i, £&),&„(£, ek)] in 
terms of the initial condition [a n (t, ek = 0),6„(t, ek = 0)]. The group transforma- 
tion corresponding to the group parameter ek can usually be written explicitly 
only for the lowest values of k. In the case of the generalized symmetries, the 
group action is obtained in principle by solving the Cauchy problem for the non- 
linear characteristic equation starting from a generic solution of the Toda lattice. 
This often cannot be done. In all cases we can construct just a few classes of 
explicit group transformations corresponding to very specific solutions of the 
Toda lattice equation, namely the solitons, the rational solutions and the peri- 
odic solutions [230]. In all cases one can use the symmetries 14.27f) to do, for 
example, symmetry reduction and to reduce the equation under consideration 
to an ordinary differential equation, or possibly a functional one. This is done 
by looking for fixed pints of the transformation, i.e. putting a n ^ k — 0, b n ^ k = 0, 
We can extend the class of symmetries by considering nonisopectral defor- 
mations of the spectral problem ll4~TH|l [65, 70, 71, 132, 162]. Thus for the Toda 
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system we obtain 



C k +H 



bn,e k J V 0>n-l — a 



a n {b r , 



>n+l, 



, r k( a n [(2n + 3)b n+ x - (2n - l)b n 
\bl - 4 + 2[(n + l)a„ - (n - lK-ii 



(4.29) 



In correspondence with ea. H4.29|) we have the evolution of the reflection coeffi- 
cient H4.24L given by 

^= M A fc+1 ii?, X €k =^X k , (4.30) 

The proof that the eqs. I|4.29(l are symmetries is done by showing that the flows 
PUf and in the space of the reflection coefficients commute. 

In addition to the above two hierarchies of symmetries (|4.27|l and (|4.29(l , we 
can construct two further cases, which, however, do not satisfy the asymptotic 
boundary conditions (|4.23l) . They are: 

£)-(!)■ fc;i)-("t':,"-t' , ) + (t)- < 4 - 31 ' 

As these exceptional symmetries do not satisfy the asymptotic boundary condi- 
tions (|4.23l) . we cannot write a corresponding evolution equation for the reflec- 
tion coefficient (|4.24|l . 

Let us now write down the lowest order symmetries for the Toda system 
H4.16JI that one can get from the hierarchies (|4.27l I4.29() . In the case of the 
Toda lattice l|2.3|) the symmetries are obtained from those of the Toda system 
(|4.16|l by using the transformation l|4.17|) . The symmetries of the Toda lattice 
and the Toda system, corresponding to the isospectral and nonisospectral flows, 
will have the same evolution of the reflection coefficient. The transformation 
(|4.17|l involves an integration (to obtain u n ). The integration constant must be 
chosen so as to satisfy the following boundary conditions: 

lim v n = 0. (4.32) 

n| — >oc 

In the case of the exceptional symmetries such an integration will provide an 
additional symmetry. 

Taking k = 0, 1, and 2 in eq. (|4.27|) we obtain the first three isospectral 
symmetries for the Toda system, namely: 

fln,eo = Q"n(b n — &n+l), 

b n ,e = a n-l ~ ««, (4.33) 



an,ei — an(bn — bf l+1 + a n _i — a„+i), 
b n ,e 1 = a n -i(b n + b n -i) - a n (b n +i + b n ), 



(4.34) 
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a«, £ 2 — a n(bn — b n+1 + a n 6 n — 2a„ + i& Jl+ i + a n -\b n -i + 2a n -\b n 

—a>n+ibn+2 — a n b n +i — 2b n + 2b n +i), 
b n ,e 2 = a n -i(b1 + b 2 l _ 1 + b n b n ~i + a„_i + a«-2 - 2) 

-a n {b 2 n + b 2 n+l + b n b n+ i + a n+1 + a n - 2). (4.35) 

The lowest nonisospectral symmetry is obtained from eq. I|4.29[l . taking k = 0. 
It is: 

an,v = a n {t(b 2 n - b 2 +1 + a n -i - a n+1 ) + (2n + 3)b n+1 - (2n - 1)6„}(4.36) 
b n . v = t{a n ^i(b n + 6 n _i) - a n {b n+ i + b n )} + b 2 n -A + 2[(n + l)a n - (n - l)a„_i]. 

The nonisospectral symmetries for k > are nonlocal. 
The exceptional symmetries <|4.31(1 are: 

a n , w = 0, b n ,[i = 1, (4.37) 
a n ,m = 2a n + ta n , b n ^ t =b n + tb n . (4.38) 

The corresponding symmetries for the Toda lattice are 

V„,e = v n (4.39) 

v n>ei =v 2 n + e^-i-"- + e »»-»»+i - 2 (4.40) 



v n , €2 =v 3 n - 2v n + e v -^ v - (*„_! + 2v„) + e 1 '"-^ 1 + 2« n ) (4.41) 

v n , v = + e - 1 -"- + e""-^ 1 - 2} - (2n - l)v n + w n (t), (4.42) 
where w n (t) is defined by the following compatible system of equations: 

w n+1 (t) - w n (t) = -2v n+1 , w n (t) = 2(e v "- Vn +i - 1). (4.43) 

Under the assumption l|4.32|l we can integrate eqs. (|4.43|) and obtain a formal 
solution. That is, we can write w n (t) in the form of an infinite sum; 

oo 

w n (t) = 2 Vj+a, (4.44) 

j=n+l 

where a is an arbitrary integration constant which can be interpreted as an 
additional symmetry. The exceptional symmetries read: 

Ufi,f*i = tv n - 2n (4.45) 



«n,M0 = t (4.46) 
and the additional one, due to the integration 

Vm-j = L ( 4 - 47 ) 



33 



4.2.1 The symmetry algebra for the Toda lattice 

To define the structure of the symmetry algebra for the Toda lattice we need 
to compute the commutation relations between the symmetries, i.e. between 
the flows commuting with the equations of the hierarchy Using the one-to-one 
correspondence between the integrable equations and the evolution equations for 
the reflection coefficients, we calculate the commutation relations between the 
symmetries and thus analyze the structure of the obtained infinite dimensional 
Lie algebra. 



If we define 



' Ak) r (ky 

i -21 *-22 



we can write the generators for the isospectral symmetries as 

= {C$[a n (b„ - 6„+i)] + C$(a n -i - a n )}d an 

+{4?M&n - b n+i)\ + ^'K-i - a n )}d bn . (4.49) 

The superscript on the generator X, is there to indicate that this is the symme- 
try generator for the Toda system (|4.1(il) . We will indicate by the superscript 
TL the symmetry generators for the Toda Lattice l|2.3|) . To these generators 
we can associate symmetry generators in the space of the reflection coefficients. 
These generators are written as 

Xl = n\ k Rd R . (4.50) 

In agreement with Lie theory, whenever R is an analytic function of e, the 
corresponding flows are given by solving the equations 

^ = ^X k R, ^- = 0, R(e k = 0) = R, \(e k = 0) = A. (4.51) 
ae k dek 

One can prove that the isospectral symmetry generators 14.4911 commute 
amongst each other 

=0, (4.52) 

by computing the corresponding commutation relation in the space of the re- 
flection coefficients 

[4 T , X%] = [fi\ k Rd R , n\ m Rd R ] = 0. (4.53) 

So far, the use of the vector fields in the reflection coefficient space has 
just reexpressed a known result, namely that the commutation of the reflection 
coefficients is rewritten as eq. Q4.53t ■ We now extend the use of vector fields in 
the space of the spectral data to the case of the nonisospectral symmetries I 
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Using the definition l|4.48|l we can introduce the generators of the nonisospectral 
symmetries for the Toda system. The symmetry vector fields are: 

Y k T = {t[C[\ +1) [a n (b n - b n+1 )} + 4* +1) K-i ~ a™)] 
+4? M(2n + 3)6„+i - (2n - l)b n )} 
+C ( $[bl - 4 + 2(n + l)a„ - 2(n - l)a„_i]R„ 
+{t[£^ +1) [a„(6„ - 6„+i)] + £^ +1) ( a ™-i ~ a «)] 

+4? [a„((2n + 3)6„+i - (2n - l)b„)] 
+4? P# - 4 + 2(n + l)a„ - 2(n - l)a n _i]H„ . (4.54) 

Taking into account eq. i|4.30[) . we can define the symmetry generators (|4.54|l 
in the space of the spectral data, as 

yl = fi\ k+ HRd R + fi 2 \ k d x . (4.55) 
Commuting y k with we have: 

[yL yl] =(m-k) [yl +m+1 - 4j>f +m _ 1 ] . (4.56) 

From the isomorphism between the spectral space and the space of the solutions, 
we conclude that the vector fields representing the symmetries of the studied 
evolution equations, satisfy the same commutation relations 

[Y fe T , Yl\ = (m - k)[Y k T +m+1 - 4yf +m _ 1 ], (4.57) 

In a similar manner we can work out the commutation relations between the 
and X m symmetry generators. We get: 

K T , S>Z] = -(1 + + MXk+m-x, (4-58) 

and consequently 

[Jtf, Y%\ = -(1 + k)Xl +m+1 + 4kX% +m _ v (4.59) 

Relations like l|4.57|l and (|4.59|) can also be checked directly, but the use of 
vector field in the reflection coefficient space is much more efficient. 

Let us now consider the commutation relations involving the exceptional 
symmetries (|4.31f) . We write them as: 

Zq = d bn (4.60) 

Zf = [2a n + ta n ]d an + [b n + tb n ]d bn . (4.61) 

As mentioned in Section 14.11 these symmetries do not satisfy the asymptotic 
conditions (|4.23|l . Hence we cannot write down the commutation relations in 
all generality for all symmetries simultaneously. We calculate explicitly the 
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commutation relations involving just Zq and Z~[ , Xq , Xj and Yq . The non 
zero commutation relations are: 

[Xq, Zj] = ~Xq, [Zq , Zj] = Zq 

K) T ) Z I\ = , {Yq T , Z\ \ = -Yq T - 8Zq, [Xj, Zq] = -2Xq, 

[X?, Zl\ = -2±l, [X T q,Yq T \ = -If, [±1,Yq T \ = -2Xl + 4l T - (4-62) 

In the case of the Toda lattice <|2.3|) we have (see eas. H4.46l I4.45|0 

Zl L = td Vn (4.63) 

Zj L = [tv n - 2n]d Vn (4.64) 

and 

Z T J( = d Vn (4.65) 

in correspondence with eq. (|4.47|) . As eqs. (|2.3|) and 14.16|l are just two different 
representations of the same system, the symmetry generators in the space of the 
spectral data are the same. Consequently the commutation relations between 
X% L and Y£ L are given by eqs. (|433l 14371 and l4~5^|l . The symmetries Xq L , 
Xf L and Yq TL , according to eas- flO^ . (H30J), (|4~4H) are given by: 

Xl h = v n d Vn , Xl L = [vl + e v ^~ v " + e v "~ v ^ - 2}d Vn 

Y TL = {t[vl t + e v ^- v " + e v "~ v ^ - 2} - (2n - l)v n , t + w n {t)}d Vn 

w n+1 (t) - w n (t) = -2v n+1 , w n {t) = 2(e Vn - Vn + 1 - 1). (4.66) 

The nonzero commutation relations are: 

\vTL ?TLi _ £TL fvTL $-TLi _ v-TL 

[Xq TL ,Yq TL ]=-XI l +u J Z tl 1 , 

[X? L , Zl L \ = -2Xl L , [Xl L ,Zl L \ = -2±l L - AZ T [ , 

[Jf[ L , Yq TL ] = -2±l L + 4X Ti + aZ^i , 

[Yq TL , Z T _\\ = pZ T \ , [Yq TL , Z T q L \ = -2Z T X L + 1 Z T _[, 

[Yq TL , Z T X L \ = ~Yq TL - 8Z^ L + 5Z T {, [Z T q L ,Z t x l \ = Z T q L , (4.67) 

where (/3, 7, S, u), a) are integration constants. The presence of these integra- 
tion constants indicates that the symmetry algebra of the Toda equation is not 
completely specified. The constants appear whenever the symmetry Y TL is 
involved. The ambiguity is related to the ambiguity in the definition of Y TL 
itself, i.e. in the solution of eq. I|4.66|l for w n (t). We fix these coefficients by 
requiring that one obtains the correct continuous limit, i.e. in the asymptotic 
limit, when h goes to zero, a combination of the generators of the Toda Lattice 
l|2.3[) and Toda system H4.16|l goes over to the symmetry algebra of the potential 
Korteweg - de Vries equation (see below). 

The commutation relations obtained above determine the structure of the 
infinite dimensional Lie symmetry algebras. The first symmetry generators are 
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given in eas. (|4.49|l . 14. 5411 . I|4.60[) . I|4.61l) and the corresponding commutation 
relations are given by eas. l|4.57l) . (|4.59() . (|4.62|) . As one can see, the symmetry 
operators Y£ and Z^ are linear in t and the coefficient of t is an isospectral 
symmetry operator X^ _ Consequently, as the operators X]: commute amongst 
each other, the commutator of with any of the Y" fc T or Z^ symmetries will 
not have any explicit time dependence and thus can be written in terms of X^ 
only. Thus the structure of the Lie algebra for the Toda system can be written 
as: 

L = L 5Li, L = {/ l ,e,/,Y 1 T ,f 2 r ,---},L 1 = {X T ,l 1 T ,...} (4.68) 

where {h — Z±,e — Zj , / = Y$ + 4Z(f} denotes a sl(2,R) subalgebra with 
[h, e] — e, [h, f] = — /, [e, /] = 2h. The algebra L is perfect, i.e. we have 
[Lq.Lq] = Lq. It is worthwhile to notice that ZQ,Zf and Xq are point sym- 
metries while all the others are generalized symmetries. 

For the Toda lattice equation the point transformations are Xq L , Zq L and 
Zi L , as for the Toda system, plus the additional Z?_\. Taking into account eqs. 
l|4.63H4.67|l . the structure of the Lie algebra is the same as that of the Toda sys- 
tern with L = {Z T [ , Zj L , Zf L , F Ti , Y 1 Ti , Y™, ■■■},L 1 = X? L , ■ 

4.2.2 Contraction of the symmetry algebras in the continuous limit 

It is well known [27, 30, 136, 137, 230], that the Toda equation has the potential 
Korteweg-de Vries equation as one of its possible continuous limits. In fact, by 
setting 

v n (t) = ~^hu(x, t) x = (n - t)h t = --^h 3 t (4.69) 
we can write eq. (|2.3|l as 

K - <W - 3u 2 x ) x = 0(h 2 ) (4.70) 

i.e. the once differentiated potential Korteweg-de Vries equation. Let us now 
rewrite the symmetry generators in the new coordinate system defined by l|4.69|l 
and develop them for small h in Taylor series. We have: 

X? L = {-u x (x,T)h- ^u T (x,T)h 3 }d u (4.71) 

X( L = {-2 Mx (x,r)/ l ~i UT (x,r)/ l 3 +0{h 5 )}d u (4.72) 

Xl L = {-Au^x^h-^Urix.r)^ +0{h 5 )}d u (4.73) 
f Ti = {2[u(x,T) + xu x (x,T) + 3TU T (x,T)] + 0(h)}d u (4.74) 
= ~d u , Z^ = ^rd u (4.75) 

Zl L = {-^T+^[x + 6Tu x (x,T)]+0(l)}d u . (4.76) 
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To obtain eqs. H4.72tBTH) we have used the following evolution for u: 



u T = u xxx + 3u 2 x . (4.77) 

The point symmetry generators written in the evolutionary form, for the poten- 
tial Korteweg-de Vries equation (|4.77l) read: 

P = u T d u , Pi=u x d Ul B = [x + 6ru x ]d u , (4.78) 
D = [u + xu x + 'Sru T ]d u , f = d u , (4.79) 

and their commutation table is: 





Po 


Pi 


B 


b 


f 


Po 








-6Pi 


-3P 





Pi 







-f 


-A 





B 









IB 





b 











-f 


f 














(4.80) 



Let us now consider the continuous limit h — > of the symmetry algebra 
of the Toda equation. We can write the simplest symmetry generators of the 
Toda equation as a linear combination of the generators l|4.71hl4~7SI) , so that in 
the continuous limit they go over to the generators of the point symmetries of 
the potential Korteweg-de Vries equation (|4.78l 14.791) : 

= ^(2l Ti -lP), P 1 = -ix Ti , D = ±Y TL , (4.81) 

B = ^(2^ + ^), f = -^f. (4.82) 

Taking into account the commutation table between the generators Xq L , Xi L , 
Z^i, Zq L Z^ L and Y Ti , given by 14.67|l and the continuous limit of Xj L 
given by eq. (|4.73|l . we get: 





Po 


Pi 


B 




b 


f 


Po 








-6Pi + 0(h 2 ) 


-3P 


+ 0(h 2 ) 





Pi 







-r + o(h 2 ) 


-Pi 


±<D(h 2 ) 





B 









2B H 


-0{h 2 ) 





b 













-f 


f 
















(4.83) 



The results contained in table l|4.83|l are obtained by setting (3 = —2, 2j + 5 = 
0, and lu = a = 0. Thus, we have reobtained in the continuous limit, all 
point symmetries of the potential KdV equation. The limit partially fixes the 
previously undetermined constants in eq. I|4.67(l . To get all point symmetries 
of the potential KdV equation we needed not only the point symmetries Xq L , 
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Zj L , Zzi and Zf L of the Toda equation, but also the higher symmetries X^ L , 

vTL 

z o ■ 

This procedure can be viewed as a new application of the concept of Lie al- 
gebra contractions. Lie algebra contractions were first introduced by Inonu and 
Wigner [105] in order to relate the group theoretical foundations of relativistic 
and nonrelativistic physics. The speed of light c was introduced as a parameter 
into the commutation relations of the Lorentz group. For c — > oo the Lorentz 
group "contracted" to the Galilei group. Lie algebra contractions thus relate dif- 
ferent Lie algebras of the same dimension, but of different isomorphism classes. 
A systematic study of contractions, relating large families of nonisomorphic Lie 
algebras of the same dimension, based on Lie algebra grading, was initiated by 
Moody and Patera [180]. 

In general Lie algebra and Lie group contractions are extremely useful when 
describing the mathematical relation between different theories. The contrac- 
tion parameter can be the Planck constant, when relating quantum systems to 
classical ones. It can be the curvature A: of a space of constant curvature, which 
for k — ► goes to a flat space. The contraction will then relate special functions 
defined e.g. on spheres, to those defined in a Euclidean space [110]. 

In our case the contraction parameter is the lattice spacing h. Some novel 
features appear. First of all, we are contracting an infinite dimensional Lie 
algebra of generalized symmetries, that of the Toda lattice. The contraction 
leads to an infinite dimensional Lie algebra, not isomorphic to the first one. 
This " target algebra" is the Lie algebra of point and generalized symmetries of 
the potential KdV equation. A particularly interesting feature is that the five 
dimensional Lie algebra of point symmetries of the potential KdV is obtained 
from a subset of point and generalized symmetries of the Toda equation. This 
5 dimensional subset is not an algebra (it is not closed under commutations). 
It does contract into a Lie algebra in the continuous limit. 

4.2.3 Backlund transformations for the Toda equation 

In addition to symmetry transformations presented in Section l|4.2.1(l . the Toda 
system admits Backlund transformations [25,29,33,49,91,122,126,133]. They 
are discrete transformations (i.e. mappings) that starting from a solution, pro- 
duce a new solution. Backlund transformations commute amongst each other, 
allowing the definition of a soliton superposition formula that endows the evo- 
lution equation with an integrability feature. Using the spectral transform [33] 
we can write down families of Backlund transformations. They are obtained by 
requiring the existence of two essentially different solutions to the Lax equa- 
tions Q4.2I I4.3JI . ip(x,t; A) and ip(x,t; A). These two solutions will be associated 
to two different solutions to the nonlinear partial differential equation l|4.1|) . 
u(x,t) and u(x,t) and consequently two different Lax pairs (L(u), M(u)) and 
(L(w), M(u)). If a transformation exists between u and u than there must exist 
a transformation between ip(x,t]X) and ip(x,t;X) and between [L{u), M(u)) 
and (L(u), M(u)). This implies that there will exist an operator D(u,u), often 
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called the Darboux operator which will relate ip(x, t; A) and ip( x i % ty> i- c - 

ip(x, t; A) = D(u, u)ip(x, t; A). (4.84) 

Taking into account the Lax equations for ip(x, t; A) and ip(x, i; A), we get from 
<|4.84[1 the following operator equations for D: 



L(u)D(u,u) = D(u,u)L(u), 

D t (u, u) = D(u, u)M{u) - M{u)D{u, u). 



(4.85) 
(4.86) 

From eqs. (14.851 14.86fl we get a class of Backlund transformations, which we 
will symbolically write as Bj(u(x,t),u(x,t), . . .) characterized by a recursion 
operator A. 

In the discrete case, the technique is basically the same, and in the case of 
the matrix discrete Schrodinger spectral problem, a generalization of the scalar 
problem (|4.18|) , the appropriate developments are found in [29] . Specializing to 
the scalar case, the class of Backlund transformations associated to the Toda 
system (|4.16|) is given by 

fa(n) - a(n)\ _ ( fl^Tr 1 {n+l)(b{n) - 6(n+l)) \ , , 
7(A) \b{n) - b(n) ) ~ 6{A > {fi(n-l)Il-\n) - fi^-^n+l) J (4 ' 87) 

where j(z) and S(z) are entire functions of their argument and we have denoted 

oo oo 

II(n) = l[a(j), fl(n) = jj o(j). (4.88) 

j=n j=n 

Above, A is the recursion operator 

p(n)b(n+l) + a(n)[q(n) + q{n+l)\ + T,(n)[b(n) - b(n+l)} 

OO 

+[<z(n) - a(n)]J2p{j) 



A 



p(n) 
q(n) 



p(n) + b(n)q(n) - E(n) + E(n-l) + [b(n) - b(n)] VVj) 



j=n 



and 



E(n) = n(n) 



(4.89) 



(4.90) 



In [30] it is proven that whenever (a n ,b n ) and (a n ,b n ) satisfy the asymptotic 
conditions (|4.23|) and the Backlund transformations l|4.87|l . the reflection coef- 
ficient satisfies the equation 



j(X)-S(X)/z 



(4.91) 
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For the Toda lattice 1)2. 3|l the one-soliton Backlund transformation, when 
the functions 7(A) and (5(A) are constant, reads: 

u{n) - u(n+l) = /3{e M ( n+1 )-«("+ 1 ) - e »(n-i)- u (n) - e "(")-"( Tl ) + e «(n)-«(«+i) j. 

(4.92) 

Equation l|4.92[l can be interpreted as a (nonlinear) three point difference equa- 
tion for u, where u is some chosen solution of (12. 311 . 

The formulas l|4.87tFOJT|) provide also much more general transformations, 
i.e. higher order Backlund transformations. If the arbitrary functions 7(A) 
and <5(A) are polynomials, then we have a finite order Backlund transforma- 
tion that can be interpreted as a composition of a finite number of one-soliton 
transformations. In more general cases, when 7(A) and (5(A) are entire functions, 
we face an infinite-order Backlund transformations. 

In the following section we discuss how Backlund transformations are related 
to continuous symmetry transformations, allowing, albeit formally, an integra- 
tion of the latter. 

4.2.4 Relation between Backlund transformations and higher sym- 
metries 

A general isospectral higher symmetry of the Toda equation is given by 

m , a n (b n - b n+1 )\ ^ (4 Q3) 



with the spectum evolution 

^=MW, f ). (4.94) 

These equations generalize eq. (|4.27|l . used above. In eqs. (|4.93l I4.94|) the 

function <j) is an entire function of its argument. Eq. (|4.94|) can be formally 
integrated in the spectral parameter (A) space, giving: 

R{\, e) = exp(/!0(A))i?(A, 0). (4.95) 

Taking into account the following definitions of A and (i in terms of z 

A = - + z, fi = - - z, /i 2 = A 2 - 4, (4.96) 

\-^,'i-±±i!, (4.97, 

we can rewrite the general Backlund transformation (|4.91l) as 

- 2-(A- M )/3(A) «5(A) 
2- (A + /x)/3(A) 7(A) 

In order to identify a general symmetry transformation with a Backlund trans- 
formation, and vice versa, we equate i?(A, e) = R(X) 
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and find that </>(A) in eq. 1)4. 95|) is given by 



1 2 - (A - (j)pW 
= ^ 2 - (A + n)(3{\) 



(4.99) 



The rhs of eq. (|4.99|) must not depend on /i. Relations l|4.96|l allow us to 
separate the exponential in (|4.98(l into two entire components E (\) and E 1 (X) 
satisfying 



exp(//0(A)) = cosh(^</)(A)) + /i 



sinh(//0(A)) 
A* 



Bo(A)+M-Si(A). (4.100) 



Noticing that the rhs of eq. (|4.99ll is rational in /i, developing /i 2 and identifying 
powers (0th and 1st) of /i, we get a system of two compatible equations 

-(2 - A/3)£ + (A 2 - 4)/3£i = -(2 - A/3), (4.101) 
-/3£ + (2 - Xf3)E 1 = (3. (4.102) 

Eqs. I|4.101l 14. 1 02|) provide us with explicit formulas relating a given general 
higher symmetry (characterized by <$>, and thus E , Ei) with a general Backhand 
transformation (characterized by 7 and <5, and thus by /3): 



/3(A) = 



2£i 



2 sinh(/i0) / ^ 



M = 

7(A) £; + A£i + l cosh(/^) + A sinh(^0)/V + 1 ' 



(4.103) 



From this equation we see that whatsoever be the symmetry, we find a Backlund 
transformation, i.e. for an arbitrary function <f> we obtain the two entire func- 
tions 7 and 5. Vice versa, given a general Backlund transformation, we can find 
a corresponding generalized symmetry 



Eq — 



2(/3 2 



1) 



A/3(2 - A/3) 



2(/3 2 - A/3 + 1) 



£1 = 



(A/3 - 2)/3 



or more explicitly, 



MX) ~ — sinh 
A* 



2(/3 2 - A/3+1)' 
(A/3 - 2)/3 



"2(/3 2 - A/3 + 1) 
In the case of a one-soliton Backlund transformation with /3 



(4.104) 

(4.105) 
1, we have: 



and we can write </>(A) as 



MX) 



X 1 

2' El = 2" 



sinh" X (VA 2 -4/2) 



Va 2 



(4.106) 



In this simple case we can write the symmetry in closed form as an infinite 
sequence of elementary symmetry transformations: 



fc=0 



(2fc)!7T 



fc!(fc-l)!2 4fc + 2 



1 k\{k+\)\ 



2 (2fc+2)! 



(4.107) 
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In this way, the existence of a one-soliton transformation implies the existence 
of an infinite-order generalized symmetry. 

Let us consider the time shift symmetry given by 0(A) = 1. Then eq. I|4.100|l 
implies that Eq = cosh// and E\ = sinh/x//x. According to (|4.103[1 the corre- 
sponding higher Backlund transformation is 

<5(A) = 2sinh fi/fj. (4.108) 
7(A) = cosh/i + Asinh/i/ju + 1 (4.109) 

This Backlund transformation, corresponding to the point symmetry studied, 
is of infinite order. 



5 Point Symmetries Transforming Solutions and 
Lattices 

5.1 Lie Point Symmetries of Ordinary Difference Schemes 
5.1.1 Ordinary Difference Schemes 

In this section we take a complementary view to the previous ones. We shall 
restrict ourselves to continuous point symmetries only. On the other hand we 
shall consider flexible difference schemes, rather than just difference equations 
on fixed lattices. For brevity and simplicity of exposition, we limit ourselves 
to ordinary difference schemes (OAS). For partial difference schemes, involving 
more than one independent variable, see Ref. [55, 143]. 

An OAS involves two objects, a difference equation and a lattice. We shall 
specify an OAS by a system of two equations, both involving two continuous 
variables x and u(x), evaluated at a discrete set of points {x n }. 

Thus, a difference scheme of order K will have the form 



WS.f.KO = 0, a = 1,2 
K = N-M + 1, n,M,NeZ, N > M, u k = u(x k ) 



(5.1) 



At this stage we are not imposing any boundary conditions, so the reference 
point x n can be arbitrarily shifted to the left, or to the right. The order K of the 
system is the number of points involved in the scheme l|5.1|) and it is assumed 
to be finite. We also assume that if the values of x k and u k are specified at 
(N — M) neighbouring point, we can calculate their values at the point to the 
right, or to the left of the given set, using equations l|5.1fl i.e., the Jacobians 
satisfy 

d(El > E2) *0, J {E ^ ] ^0. (5.2) 



d(x n +N, u h+ n) d{x n+ M,u n+ M) 

A continuous limit, when the spacings between all neighbouring points go to 
zero, if it exists, will take one of the equations (|5.1|l into a differential equation 
of order K 1 < K, the other into an identity (like = 0). 
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When taking the continuous limit it is convenient to introduce different quan- 
tities, namely differences between neighbouring points and discrete derivatives 
like 



Hs^-iXfi) X n -\-i ^ni h—{Xyi) X n X n - 



1- 



u n +i —u n u n - u n -x . . 

u x = , u SL = — , (5.3) 



V<xx — 2 , . . . 

•**n+l *£n — 1 

In the continuous limit, we have 

h + — * 0, h- — > 0, u, x — > u', Ux — > "ra — > u". 

As a clarifying example of the meaning of the difference scheme (|5.1fl . let 
us consider a three-point scheme that will approximate a second-order linear 
difference equation: 

U n+ 1 - 2u n + Un-i 

El = — -2 u n = 0, (5.4) 

[Xn+l Xn ) 

E 2 = x n+ i - 2x n + x n -i = 0. (5.5) 



The solution of eq. E 2 — 0, determines a uniform lattice 

x n = hn + xq. (5-6) 

The scale h and the origin xq in eq. (|5.6[) are not fixed by eq. 15.5fl . instead 
they appear as integration constants, i.e., they are a priori arbitrary. Once 
they are chosen, eq. (|5.4() reduces to a linear difference equation with constant 
coefficients, since we have x n +i — x n — h. Thus, a solution of eq. (|5.4|l will have 
the form 

u n = X Xn . (5.7) 

Substituting 15.711 into 15. 4|) we obtain the general solution of the difference 
scheme (|5.4ll . (15.51) . 

u(x n ) = ciX\ n + C2A2" , x n = hn + xo, (5-8) 

/ 2 + h 2 ± ftVlTT? 

Ai,2 - ^ 2 

The solution (|5.8(l of system (|5.4|) - H5.5|) depends on 4 arbitrary constants ci, C2, 
ft. and a;o- 

Now let us consider a general three-point scheme of the form 

E a (x n -i,x n ,x n+ i,u n -i,Un,'U n+ i) = 0, a =1,2, (5.9) 

satisfying (|5.2|l with iV = 1, M = — 1 (possibly after an up or down shifting). 
The two conditions on the Jacobians 13. 2H are sufficient to allow us to calculate 
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(x n+1 ,u n+1 ) if (x n -i,Un-i,x n ,Un) are known. Similarly, {x n -i,u n -i) can be 
calculated if (x n , u n , x n +i, u n+ \) arc known. The general solution of the scheme 
H5.9fl will hence depend on 4 arbitrary constants and will have the form 

U n = f{x n ,Cl,C2,Cz,Ci) (5.10) 
X n = <f>(n,Cl,C 2 ,C 3 ,Ci). (5-11) 

A more standard approach to difference equations would be to consider a 
fixed equally spaced lattice, e.g., with spacing h = 1. We can then identify the 
continuous variable x, sampled at discrete points x n , with the discrete variable 
n, 

x n = n. (5.12) 

Instead of a difference scheme we then have a difference equation, Ike the one 
we considered in Sectional 

E({u k y^ +M ) = 0, (5.13) 

involving K = N — M + 1 points. Its general solution has the form 

u n = f(n, ci,c 2 , . . .c N - M ) , (5-14) 

i.e., it depends on N — M constants. 

Below, when studying point symmetries of discrete equations we will see 
the advantage of considering difference systems like system l|5.1|l . A treasury 
of information on difference equations and their solutions can be found in the 
classical book by Milne-Thompson [179]. 



5.1.2 Point Symmetries of Ordinary Difference Schemes 

In this section we shall follow rather closely the article [141]. We shall define 
the symmetry group of an ordinary difference scheme in the same manner as 
for ODEs: it is, a group of continuous local point transformations of the form 
l|1.2|l taking solutions of the OAS l|5.1[) into solutions of the same scheme. The 
transformations considered being continuous, we will adopt an infinitesimal ap- 
proach, as in eq. (|1.3|l . We drop the labels i and a, since we are considering the 
case of one independent and one dependent variable only. 

As in the case of differential equations, our basic tool will be vector fields of 
the form 1|1.4|) . In the case of OAS they will have the form 

X = £{x,u)d x + <j>(x,u)d u , (5.15) 

with 

X = X n , V> — Un — u(Xn). 

Because we are considering point transformations, £ and (f> in (|5.15|) depend on 
x and u at one point only. 
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The prolongation of the vector field X is different from that of the case of 
ODEs. Instead of prolonging to derivatives, we prolong to all points of the 
lattice figuring in scheme l|5.1|) . Thus we set 

n+N n+N 

prX = ^2 £,( x k,u k )d Xk + ^2 </)(x k ,u k )d Uk . (5.16) 

k=n+M k=n+M 

In these terms the requirement that the transformed function u(x) should satisfy 
the same OAS as the original u(x) is expressed by the requirement 

prXE a | Bl=£2= o=0, a =1,2. (5.17) 

Since we must respect both the difference equation and the lattice, we have 
two conditions l|5.17|l from which to determine £,{x,u) and 4>{x,u). Since each 
of these functions depends on a single point (x,u) and the prolongation l|5.16|l 
introduces N — M + 1 points in space X x U, equation (|5.17() will imply a 
system of determining equations for £ and (/). Moreover, in general this will be 
an overdetermined system of linear functional equations that we transform into 
an overdetermined system of linear differential equations [9, 10]. 

Let us first of all check that the prolongations (|5.16|l have the continuous 
limit. We consider the first prolongation in the discrete case 

pr x X = ^(x,u)d x +<j)(x 1 u)d u + £(x+, u + )d x+ + <f>(x+, u + )d u+ , (5.18) 

X + = X n +1, U + = U n +X, 

and apply it to a function of the variables x, u, h + , u x (15.31) . The continuous 
limit is recovered by taking h — > and putting 

u + = u{x+) = u(x) + hu (x) + . . . , (5.19) 
£(x + ,u + ) = £(x,u) +£ <x h + £, u u'h + . . . , 
<t>(x + ,u+) — <t>(x,u) + (j> tSC h + (f> lU u'h + 

where u' is the (continuous) derivative of u(x). We have 

pr 1 D XF(x, u, h, u x ) = {£d x + <f>d u + [£, x + £, u u']hd h + 
+ [-v>'(£,x + (,uU r ) + 4>. x + 0,«w']a„ x }F. 

In the limit when h — * we obtain 

lim prjjX = £(x, u)d x + 4>{x, u)d u + (f> x (x, u, u')d u r, (5.20) 

4> x - + (<P, U - e,x)«' - iuu' 2 , (5.2i) 

which is the correct expression for the first prolongation in the continuous case 
[194]. The proof that the n-th prolongation l|5.16|l has the correct continuous 
limit can be performed by induction. 
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5.2 Examples 

1. Power nonlinearity on a uniform lattice. 

Let us consider a difference scheme that is a discretization of the ODE 

u" - u N = 0, N^0,1. (5.22) 

For N ^ —3, eq. I|5.22f) is invariant under a two-dimensional Lie group whose 
Lie algebra is given by 

X x = d x , X 2 = {N- l)xd x - 2ud u (5.23) 

(translations and dilations). For N — — 3 the symmetry algebra is three- 
dimensional, isomorphic to sl(2,R), i.e., it contains a third element in addition 
to (I5.23|) . A convenient basis for the symmetry algebra of the equation 

u" - U - 3 = (5.24) 

is 

X x = d x , X 2 = 2xd x + ud u , X 3 = x(xd x + ud u ). (5.25) 
A very natural OAS that has (|5.22(l as its continuous limit is 

Ex = -, T2 u n = AT ^ 0, 1 (5.26) 

E-2 = x n+1 - 2x n + x n -x = 0. (5-27) 

Let us now apply the symmetry algorithm described in Section 15.1.21 to 
system l|5.2bll l|5.27|l . To illustrate the method, we shall present all calculations 
in detail. 

First, we choose two variables that will be substituted in eq. I|5.17|l . once the 
prolonged vector field l|5.16|l is applied to system (|5.26l) H5.27fl . namely 

x n+1 = 2x n - x n -x (5.28) 

U n +X = (X„ - X n -x) 2 Un + 2u n - U n -X 

We apply pr X of (|5.1()|) to eq. I|5.27|l and obtain 

£(x n+ x,u n +x) - 2£(x„,it„) + £(x n -i, u n -x) = 0, (5.29) 

where, x n , u n x„_i, u n _x are independent, but x n+1 , u n+ x are expressed in 
terms of these quantities, as in eq. H5.28fl . Taking this into acccount, we differ- 
entiate H5.29fl first with respect to u n -x, then with respect to u n . We obtain 
successively 

_ £,u„ +1 (Xn+l, Un+x) + 1 (^n- 1 : M n- 1 ) = (5.30) 

(N{x n -xn-x) 2 ^- 1 +2)£, Un+lUn+1 (x n+ x,u n+ i) = 0. (5.31) 
Eq. (|5.31|) is the desired one-term equation. It implies that 

£(x,u) = a(x)u + b(x). (5.32) 
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Substituting (|S~32l) into (|5~3UI) we obtain 

-a(2x n - x n -i) + o(a; n _i) = 0. (5.33) 

Differentiating with respect to x n , we obtain a = a,Q = const. Finally, we 
substitute l|5.32[l with a = ao into (|5.29|) and obtain 

a = a = 0, b{2x n - x„_i) - 2b(x n ) + b(x n -i) = 0, (5.34) 

and hence 

£ = b = b lX + b , (5.35) 

where &o and b\ are constants. To obtain the function <f>(x n ,u n ), we apply prX 
to eq. (|5.26|) and obtain 

) - 2<j>(x n ,u n ) + <f>(x n -i,u n -i) 

-{x n - Xn^fiN^Xn,^)^- 1 + 2b lU % ) = 0(5.36) 

Differentiating successively with respect to it n _i and u„ (taking (|5.28|l into 
account), we obtain 

_ 0,u„ + i( a; ri+l,Wn+l) + 4>,u n - 1 (Xn-l,Un-l) = (5.37) 

{N{x n -x n - X fu% + 2)0,„„ +1 „„ +1 = 0, (5.38) 

and hence 

4> = (j>iu + 4>o(x), 0i = const. (5.39) 
Eq. (|5.36|) now reduces to 

4> (2x n - x n -i) - 24> (x n ) + <p {x n -i) 

-{x n - x n ^) 2 ((N ~ l)0x + 26i)< 

-iV(a: n - ar„_i) 2 0o^ _1 = 0- (5.40) 

We have TV ^ 0, 1 and hence l|5.40|l implies that 

0o = O, (N — l)0i + 2bi = 0. (5.41) 

We have thus proven that the symmetry algebra of the OAS H5.26fl - I|5.27|) is 
the same as that of the ODE l|5.22|l . namely the algebra Q5.23[) . 

We observe that the value A*" = — 3 is not distinguished here and that system 
l|5.26[l - l|5.27jl is not invariant under SL(2,M) for N = -3. Actually, a difference 
scheme invariant under SX(2,R) does exist and it will have eq. I|5.24|l as its 
continuous limit. It will not however have the form H5.26fl - i|5.27|) . an( l the 
lattice will not be uniform [56,57]. The corresponding SL(2,R) invariant scheme 
is presented below in Section [5.3. 31 

Had we taken a two-point lattice, x n+ \ — x n = h with h fixed, instead of 
E2 = as in 15.27|l . we would only have obtained translational invariance for the 
equation (|5.26|) and lost the dilational invariance represented by X2 of eq. (|5.23|) . 
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2. An OAS involving an arbitrary function on a uniform lattice 
We consider 

u n+ i - 2u n + u„-i /n„-u n _i\ . . 

m = — -, — ^ / I — I = u, (,o- 4 ^j 

E 2 = x n+ i - 2x n + a;„_! = 0, (5.43) 

where f(z) is some sufficiently smooth function satisfying 

f"{z) $ 0. (5.44) 

The continuous limit of eq. 15.42fl and l|5.43|) is 

u" - f{u') = 0, (5.45) 

and it is invariant under a two-dimensional group with Lie algebra, 

X 1 = d x , X 2 =d u , (5.46) 

for any function f(u'). For certain functions / the symmetry group is three- 
dimensional, where the additional basis element of the Lie algebra is 

X A = (ax + bu)d x + (cx + du)d u . (5-47) 

Now let us consider the discrete system I|5.42[) - (I5.43[) . Before applying prX 
to this system we choose two variables to substitute in eq. (|5.17|) . namely 

x n+ i = 2x n - x n -i (5.48) 

1 i / \2 r ( U n ~ u n— 1 

Wn+1 = 2u„ — U„_i + (X n — X n -i) f I 

\ %n %n— 1 

Applying pr X to eq. (|5.43|) and performing the same kind of passages necessary 
to solve the determining equation as we did for eq. (|5.28|l we get 

£ = ax + f3, (5.49) 

with a — const, (3 = const. Now let us apply pi X to eq. (|5.42|l . We get the 
following equation for <j) 

4>(x + h,u + hz + h 2 /(z)) — 2(f>(x, u) + <p(x — h,u — hz) 

= 2ah* f(z) + f'(z) ^MJpMlM _ M ) . (5.50) 

where we have defined 

z = — — h = x n+ i-x n . (5.51) 

In general, eq. (|5.50() is quite difficult to solve and for most functions f(z) it 
has no further solutions. 
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However, if we make the choice 

f(z) = e~ z (5.52) 
we find that eq. (|5.50|l is solved by putting 

a = 1, (f) = x + u. (5.53) 
Finally, we find that the system 

_ = e .„^„_ lj Xn+1 _ 2 Xn + Xn l = o (5.54) 

is left invariant by a three-dimensional transformation group, generated by the 
solvable Lie algebra with basis 

Xi = d x \ X 2 = d u ; X 3 = xd x + (x + u)d u . (5.55) 

For further examples see the original article [141] and the lectures [238]. 

5.3 Symmetry Preserving Discretization of Ordinary Dif- 
ferential Equations 

5.3.1 General Comments 

In Section IS~T1 we assumed that an OAS (|5.1|l is given and we showed how to 
determine its symmetries. 

Here we will discuss a different problem, namely the construction of OAS 
with a priori given symmetry groups. More specifically, we start from a given 
ODE 

E(x,y,y,y,...)=0 (5.56) 

and its symmetry algebra, realized by vector fields of the form (|5.15|) . We now 
wish to construct an OAS (| f> . 1 1) . approximating the ODE H5.56(l and having the 
same symmetry algebra (and the same symmetry group). 

This can be done systematically, once the order of the ODE (|5.56|) is fixed. 
In general, the motivation for such a study is multifold. In physical applications 
the symmetry may actually be more important than the equation itself. A 
discrete scheme with the correct symmetries has a good chance of describing the 
physics correctly. This is specially true if the underlying phenomena really are 
discrete and the differential equations come from a continuum approximation. 
Furthermore, the existence of point symmetries for differential and difference 
equations makes it possible to obtain explicit analytical solutions. Finally, a 
discretization respecting point symmetries should provide improved numerical 
methods. 

Let us first outline the general method of discretization. If the ODE l|5.56[) 
is of order N we need a OAS involving at least N + 1 points, i.e. N + 1 pairs 

{a:i, W ;i = l,... > JV + l}. (5.57) 

The procedure is as follows 
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1. Take the Lie algebra g of the symmetry group 6 of the ODE (|5.56ll and 
prolong the (known) vector fields {X\, . . . , X n } to all N + 1 points (|5.57() . 
as in eq. (|5.16|) . 

2. Find a basis for all invariants of the (prolonged) Lie algebra g in the space 
l|5.57(l of independent and dependent variables. Such a basis will consist 
of K functionally independent invariants 

I a =I a (x 1 ,...,x N+1 ,yi,...,y N+1 ), l<a<K. (5.58) 

They are determined by solving the differential equations 

prX i I a (x 1 ,...,x N+1 ,y 1) ... ) y N+1 ) = 0, i=l,...,n. (5.59) 

The actual number K satisfies 

Jf = 2iV + 2-(dimg-dimg ) (5.60) 

where g is the Lie algebra of the group C stabilizing the N + 1 
points (|5.57|) . 

We need at least two independent invariants of the form (|5.58f) to write 
an invariant difference scheme. 

3. If the number of invariants is not sufficient, we can make use of invariant 
manifolds. To find them, we first write out the matrix of coefficients of 
the prolonged vector fields {Xi, . . . , X n } : 

£lJV+l 011 012 • ■ ■ 01JV+1 \ 

: : : : (5.61) 

S.nN+1 <Pnl 0n2 ■ ■ • 4>nN+l J 

and determine the manifolds on which the rank of M, Rank(M), satisfies 

Rank(M) < min(n, 2N + 2), (5.62) 

i.e. is less than maximal. The invariant manifolds are then obtained by 
requiring that eq. I|5.59|l be satisfied on the manifold satisfying eq. (|5.62|) . 

5.3.2 Symmetries of Second Order ODEs 

Let us now restrict to the case of a second order ODE 

u = F(x, u, it). (5.63) 

Sophus Lie gave a symmetry classification of second order ODE's (over the field 
of complex numbers C ) [163, 164]. A similar classification over M. is much more 
recent [173,174]. 

The main classification results can be summed up as follows. 



M = 



V?nl ?n2 
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1. The dimension n = dimg of the symmetry algebra of eq. H5.63[) can be 
dimg = 0,1,2,3 or 8. 



2. If we have dimg = 1 we can decrease the order of eq. (I5.63[) by one. If the 
dimension is dimg > 2 we can integrate by quadratures. 

3. If we have dimg = 8, then the symmetry algebra is sZ(3, C), or sZ(3, R), 
respectively. The equation can be transformed into y — by a point 
transformation. 

Further symmetry results are due to E. Noether [191] and Bessel-Hagen [19]. 
Every ODE l|5.63|l can be interpreted as an Euler-Lagrange equation for some 
Lagrangian density 

C = C(x,u,u). (5.64) 

The equation is 

^-^)=0, ^ = ^+4 + 4 + .... (5.65) 

OIL Oil OX OU OU 

An infinitesimal divergence symmetry, or a Lagrangian symmetry is a vector 
field X (|5.15|l satisfying 

prX{£) + CD(£) = D{V), V = V{x, u), (5.66) 

where V is some function of x and u. A symmetry of the Lagrangian C is 
always a symmetry of the Euler-Lagrange equation (|5.65|) . however equation 
H5.65f) may have additional, non-Lagrangian symmetries. 

A relevant symmetry result is that if we have dimg = 1 , or dimg = 2 for eq. 
H5.63J) . then there always exists a Lagrangian having the same symmetry. For 
dim£ = 3, at least a two-dimensional subalgebra of the Lagrangian symmetries 
exists. For dimg = 8 a four-dimensional solvable subalgebra of Lagrangian 
symmetries exists. 



5.3.3 Symmetries of the Three-point Difference Schemes 

A symmetry classification of three-point difference schemes was performed quite 
recently [56,57]. It is similar to Lie's classification of second order ODE's and 
goes over into this classification in the continuous limit. We shall now review 
the main results of the classification following the method outlined in Section 
l5.3.1l above. 

Sophus Lie [165] gave a classification of all finite dimensional Lie algebras 
that can be realized by vector fields of the form H5.15J) . This was done over the 
field C and thus amounts to a classification of finite dimensional subalgebras 
of diff(2,C), the Lie algebra of the group of diffeomorphisms of the complex 
plane C 2 . A similar classification of finite dimensional subalgebras of diff(2, R) 
exists [81], but we restrict ourselves to the simpler complex case. 

We use the following notation for 3 neighbouring points on the lattice: 

= Xn-i, X = X n , X + = X n+ i, U_ = Un-i, u = u n , u + = u n+1 . (5.67) 
Let us now proceed by dimension of the symmetry algebras. 
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dimg =1 A single vector field can always be rectified into the form 

Q 

A 1A : X x = — (5.68) 
ou 

The invariant ODE is 

u = F(x,u). (5.69) 
Putting u = y we obtain a first order ODE. 
The difference invariants of X\ are 

x, h+ = x+ — x, h- = x — X- , rj + = u+ — u, rj^ = u — u_. (5.70) 

Using these invariants and the notation (|5.3|l . we can write a difference 
scheme 

y xx = F(x, —^=,h-), h + = h-G(x, — — (5.71) 

This scheme goes into eq. H5.69(l if we require that the otherwise arbitrary 
functions F and G satisfy 

lim F(x, x — — , h-) = Fix, u), lim G(x. x ~^ — — ,/i_)<oo. 

(5.72) 

dimg = 2 Precisely four equivalence classes of two-dimensional subalgebras of 
diff(2, C) exist. Let us consider them separately. 

-42.1 

X 1 = d x , X 2 = d u (5.73) 

The algebra A 2i i is Abelian, the elements X\ and X 2 are linearly non- 
connected (linearly independent in any point [x, y)). The invariant 
ODE is 

il = F(u), (5.74) 
and can be immediately integrated. 

An invariant difference scheme is given by any two relations between 
the invariants h + , h—, ?y+, r\- of eq. (|5.7U|I . for instance 

u xx = F(^-^,h-), h+ = h-Gi-Z-j-Z, h-), (5.75) 
with conditions l|5.72|l imposed on the functions F and G. 

A 2 ,2 

Xi=d a , X 2 = xd x + ud u (5.76) 

This Lie algebra in non- Abelian, the two elements are linearly non- 
connected. The invariant ODE is 

u=-F(u). (5.77) 
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A basis for the difference invariance is 

h + /i 

{xu xx , u x +u x , — , — } (5.78) 
~~ — h- x 

so a possible invariant difference scheme is 

1 u x + u^ h-. ,u x +u x /i_ 
u x£ =-F( , — ),h+ = h-G{ , — ). (5.79) 

X £ X £ X 

A 2 .3 

Xi = d u , X 2 = xd u (5.80) 

The algebra is Abelian, the elements X\ and X 2 are linearly con- 
nected. The invariant ODE is 

u = F{x). (5.81) 

This equation is linear and hence has an eight dimensional symmetry 
algebra (of which A 2:3 is just a subalgebra). 
The difference invariants are 

{u xx , x, h+, h-} (5.82) 

so the invariant difference scheme will also be linear (at least in the 
dependent variable u). 

A 2A 

Xi = d u , X 2 = ud u (5.83) 

The algebra is non- Abelian and isomorphic to A 2 , 2l but with linearly 
connected elements. The invariant ODE is again linear, 

u = F(x)it, (5.84) 



as is the invariant difference scheme. Eq. Ij5.84|l is invariant under 
the group 5L(3,C). Difference invariants are 

U = 2 U 2 , x, h+, h-} (5.85) 

1l x "T Uj^ 

and a possible invariant OAS is 

2 Ux ~ = F(x,hJ), G(x,h+,h-) = 0. (5.86) 
u x -\- 

dimg = 3 We now turn to difference schemes invariant under three-dimensional 
symmetry groups. We will restrict ourselves to the case when the corre- 
sponding ODE is nonlinear. Hence we will omit all algebras that contain 
t4 2j 3 or A 2: 4 subalgebras (they were considered in [56]). 
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4m 

Xi=d x , X 2 = d u , X 3 = xd x + kud u , k^O,^, 1,2 (5.87) 
The invariant ODE is 

u^u^t. (5.88) 

For fc = 1 there is no invariant second order equation; for k = 2 the 
equation is linear, for k = | it is trasformable into a linear equation 
and has a larger symmetry group. 
Difference invariants are 

h = ^, h=u x h}+ k , I 3 =u s ht k . (5.89) 
A simple invariant difference scheme is 

u xx = ( ^) l T— y/( ^fe- *), - 2^/i- )• 

(5.90) 

We shall see below that other invariant schemes may be more conve- 
nient. 

4,2 

X 1 =d x , X 2 =d u , X 3 = xd x + (x + u)d u , (5.91) 
The invariant ODE is 

u = e~". (5.92) 
Difference invariants in this case are 

h = ^, h = h+e- u °, h = h-e- u *-. (5.93) 
A possible invariant scheme is 

Mix = e t~ f (y/ h-h+e h+ = h-g{y/h-h+e ^). 

(5.94) 

No further solvable three-dimensional subalgebras of cliff (2, C) exist 
(though there is another family for diff(2,K) [56]). 

Two inequivalent realizations of si (2, C) exist. Let us consider them sep- 
arately. 

4,3 

Xi=d x , X 2 = 2xd x + ud u , X 3 = x 2 d x + xud u , (5.95) 
The corresponding invariant ODE is 

u = u~ 3 , (5.96) 
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and its general solution is 

1 



u 2 = A(x - x ) 2 + -1 A^O. (5.97) 



A convenient set of difference invariants is 



h = h = r^-- U -±, (5.98) 
^3 = t — ; ) / 4 = 



Any three of these are independent; the four satisfy the identity 

hh = hh- (5.99) 
An invariant difference scheme can be written as 
2(/ 2 +J 3 -l) = ilh-^-^Iihh), h+h = ihl 2 g(hl 2 ), (5.100) 

^3 

i.e. 

1 1 ,h+ . h-. . 1 h + h- 

u xx — t — —r — o( — + — /bt — rr - )' 5 - 101 

h + ^ h— 4 h + h_ . 1 h + h- 
U+ u- u h+ + h- u 2 h+ + h- 

For / = g = 1 this scheme approximates the ODE Q5.96[) . 

-43,4 

X\ = d u , X 2 = xd x + ud u , X 3 = x 2 d x + {-x 2 + u 2 )d u . 

(5.102) 

This algebra is again sZ(2,C) and can be transformed into 

Yi = 9 x + d u , Y 2 = xd x + ud u , Y 3 = x 2 d x + u 2 d u . (5.103) 

The realization (|5.1U3|) (and hence also 1(5.102(1 1 is imprimitive; 
(15.91(1 is primitive. Hence A3. 4 and A3. 3 are not equivalent. The 
invariant ODE for the algebra (|5.1U2|I is 

xu = C(l +ii 2 )i + ii 2 ), (5.104) 

where C is a constant. The general integral of eq. 1(5.104(1 can 
be written as 

(x-x ) 2 + (u-u ) 2 = (^) 2 , (5.105) 

x 2 + {u-u ) 2 = x 2 ,C = Q, (5.106) 
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where xq and uq are integration constants. 

The difference invariants corresponding to the algebra l|5.103|l 
are 

j 1 = £+Zf (!+„!), /2 ^^^(i + M 2 }j (510?) 



(x + — x)(x — X-) ,., 2 . 
h = - , -{{h+ul + x+ + x)u s _ + 

+ {h-u 2 x — x- — x)u x }. 

An invariant scheme representing the ODE (|5.1()4Jl can be writ- 
ten as 

h = C(!±±k)i, h = I 2 , (5.108) 
(this is not the most general such scheme). 



5.3.4 Lagrangian formalism and solutions of three-point OAS 

In Section 15.3.21 we presented a Lagrangian formalism for the integration 
of second order ODE's. Let us now adapt it to OAS [57]. 

The La grangian density (|5.(j4[) will now be a two-point function 

C = C{x, u, x+, u + ). (5.109) 

Instead of the Euler-Lagrange equation H5.(i5(l we have two quasiextremal 
equations [51,54,57] corresponding to "discrete variational derivatives" of 
C with respect to x and y independently 

sc = h dC dC~ =Q 
5u + du du 

where C~ is obtained by downshifting C (replacing n by n— 1 everywhere). 
In the continuous limit both quasiextremal equations reduce to the same 
Euler-Lagrange equation. Thus, the two quasiextremal equations together 
can be viewed as an OAS, where e.g. the difference between them defines 
the lattice. 

The Lagrangian density (|5.109|l will be divergence invariant under the 
transformation generated by vector field X, if it satisfies 

P rX(C)+CD+(0=D + (V) (5.111) 

for some function V(x,u) where D + (f) is the discrete total derivative 

D + f(x,u) = f(* + h,u(x + h))-f(x,u) (5 n2) 
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Each infinitesimal Lagrangian divergence symmetry operator X will pro- 
vide one first integral of the quasiextremal equation 

dCr dCr 
h-4>s- + h-$^— + £C- - V = K (5.113) 
ou ox 

[57]. These first integrals will have the form 

f a (x,x + ,u, u + ) = K a , a = l, .... (5.114) 

Thus, if we have two first integrals, we are left with a two point OAS to 
solve. If we have three first integrals, then the quasiextremal equations 
reduce to a single two point difference equation, e.g. involving just x n and 
x n+ \. This can often be solved explicitly [179]. 

This procedure has been systematically applied to three point OAS in the 
original article [57]. For brevity we will just consider some examples here. 

Let us first consider a two-dimensional Abelian Lie algebra and the corre- 
sponding invariant second order ODE: 

Xi = d x , X 2 = d u , u = F{u). (5.115) 

This equation is the Euler-Lagrange equation for the Lagrangian 

C = u + G(ii), G=^, (5.116) 
r 

and both symmetries are Lagrangian ones 

prXtL + CD(Zi) = 0, P rX 2 C + CD{^ 2 ) = 1 = D{x). (5.117) 

The corresponding two first integrals are 

Jx = u + G(u) - uG(u), J 2 = G(u) - x (5.118) 

Introducing H as the inverse function of G we have 

y = H(J 2 + x), H{J 2 + x) = [G}-\J 2 + x). (5.119) 

Substituting into the first equation in eq. (|f> . 1 1 H|) . we obtain the general 
solution of l|5. 115(1 as 

y{x) =Ji- G[H{J 2 + x)] + (J a + x)H{J 2 + x). (5.120) 

Now let us consider the discrete case. We introduce the discrete analogue 
of i5Tmi) as 

c= u + u± +G{Ux) (5J21) 
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for some smooth function G. Equations H5.117[) hold (with D interpreted 
as the discrete total derivative D + ). The two quasiextremal equations are 



X+ 2 X ~ - G(u x ) + G( Ug ) = 0, (5.122) 

u x G(u x ) - UxG(Mx) - G(u x ) + G?(«x) - " + 2 "~ = 

The two first integrals obtained using Noether's theorem in this case can 
be written as 

G(u x ) - ^± = B (5.123) 
-u x G(u x ) + G(u x ) + u + - x)u x = A. (5.124) 

In principle, these two integrals can be solved to obtain 

u x = H[B + ~(x + + x)], u = $(A,B,x,x+), (5.125) 

where H[z] = [G]^ 1 (z) and $ is obtained by solving eq. I|5.124|) . once 
u x = H is substituted into this equation. A three point difference equation 
for x n+ 2, x n+i an d x n , not involving u is obtained from the consistency 
condition u x = ""^1"" ■ In general this equation is difficult to solve. We 
shall follow a different procedure which is less general, but works well when 
the considered OAS has a three dimensional solvable symmetry algebra 
with {d x ,d u } as a subalgebra. We add a third equation to the system 
namely 

x + — x 

X — X- 

The general solution of eq. (|5.I26|I is 



= l + e. (5.126) 



x n = (x + B){1 + e) n - B (5.127) 

where xo and B are integration constants. We will identify B with the 
constant in eq. I|5.I23|I . but leave e as an arbitrary constant. Eq. H5.I27fl 
defines an exponential lattice (for f / 0). Using l|5.f 27(1 together with 
HI and Om we find 

u x = H[(x n +B)(1 + |)], H[z] = {G}- 1 ^) (5.128) 
u n = A + (x n + B)H[(x n + + |)] - G[H(x n + B)(l + |)] (5.129) 

There is no guarantee that equation (|5.f 28|) and (|5.f 29|) are compatible. 
However, let us consider the two special cases with three-dimensional solv- 
able symmetry algebras, namely algebras A 3: i and A 3j2 of Section IB. 3. 31 
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Algebra ^3^1. We choose G(u x ) to be 

G(u x )= (fc ~ 1} uf*, k ^ 0,1 (5.130) 
From eq. ((5.128(1 and l(5.129|) we obtain 

«, = (-^ I ) fe - 1 4- 1 (l + |) fe - 1 (5.131) 
u n = lij^j)"' 1 ^ + B)\l + ^[l + (1 - fc)|] (5.132) 

The consistency condition (for u x to be the discrete derivative of u n ) 
provides us with a transcedental equation for e: 

[(l + £) fc -l][l + (l-£;)|] =fcc (5-133) 

In the continuous limit we take e — > and w n given by eq. I|5.132|l 
goes to the general solution of the ODE 1)5.86(1 . In eq. ((5.133(1 terms 
of order e , e, and e 2 cancel. The solution u n coincides with the 
continuous limit up to terms of order e 2 . 

We mention that in the special case fc — — 1 all three symmetries of 
the OAS are Lagrangian ones and in this case eq. 1(5.1331) is identi- 
cally satisfied for any e. 

Algebra A 3 2 - We choose G(u x ) to be 

G{u x )=e u * (5.134) 

and obtain 

u x = ln(x„ + S)(l + |), (5.135) 

u n = (x n + B) ln(.T„ + B) + A + (5.136) 
+ ( Xn + B)[\n{\+ e -)-(l+ ( -)]. 

The expressions ((5.135(1 and 1(5.136(1 are consistent if e satisfies 

e(l + |)- (l + e)ln(l+e) =0 (5.137) 

Again ((5.136(1 coincides with its continuous limit up to terms of order 
e 2 and in ((5.137(1 terms of order e°, e 1 and e 2 cancel. 

For the sZ(3, K) algebra ^3,3 all three symmetry operators X\, X2 and 
X3 correspond to Lagrangian symmetries. The corresponding OAS is 
integrated in [57]. 
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6 Conclusions and Open Problems 



Let us compare the symmetry approach for difference equations with that for 
differential ones. In both cases one is interested in transformations taking solu- 
tions into solutions and in both cases, choices have to be made. Lie's choice of 
point symmetries of differential equations is so natural, that it is often forgotten 
that it is also just an Ansatz : the vector fields should depend on the indepen- 
dent and dependent variables only. Generally speaking, this Ansatz has turned 
out to be the most useful one. Contact symmetries have much fewer applica- 
tions than point ones. Generalized symmetries are mainly of use for identifying 
integrable nonlinear partial differential equations. 

In this review we have stressed the fact that for difference equations this 
choice of pure point transformations must be modified. Without significant 
modifications it remains mainly fruitful for differential-difference equations rather 
than for purely difference ones (see Section [5Jl. 

Which modifications are needed for difference equations depends on the ap- 
plication that we have in mind. For differential equations there are two main 
types of applications in physics. In the first, the equations are already known 
and group theory is used to solve them. In the second, the symmetries of the 
physical problem are known and are used in building up the theoretical model, 
i.e. the symmetries precede the equations. These two aspects are also present 
in the case of difference equations, but there are new features. First of all, the 
physical processes that are being described may be discrete and the lattices in- 
volved may be real physical objects. If we are considering linear theories, like 
quantum mechanics, or quantum field theory on a lattice, then the generalized 
point symmetries of Section [3] are extremely promising. A mathematical tool, 
umbral calculus, is ready to be used, both to solve equations and also to formu- 
late models. For nonlinear theories on given fixed lattices the most appropriate 
symmetry approach involves generalized symmetries, as reviewed in Section 
Their main application is as in the continuous case: to identify integrable sys- 
tems on lattices. Moreover it can be used to get new interesting solutions. An 
interesting feature is that some point symmetries of differential equations, in 
particular dilations, appear as generalized symmetries for difference equations. 

The second type of application of difference equations in physics is more 
practical and in a certain sense, less fundamental. We have in mind the situ- 
ation when the physical processes are really continuous and are described by 
differential equations. These are then discretized in order to solve them. The 
lattices used are then our choice and they can be chosen in a symmetry adapted 
way. Moreover, as shown in Section [5] the difference equation and lattice are 
both part of a "difference scheme" and the actual lattice is part of the solution 
of this scheme. We can then restrict ourselves to point transformations, but 
they act simultaneously on the solutions and on the lattice. 

In an attempt to keep this review reasonably short, we have left out many 
interesting and important topics. Among them we have not included a complete 
discussion of partial difference equations on transforming symmetry-adapted 
lattices [55,58, 143] and the use of Lie point symmetries to get conditions for the 
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linearizability of difference equations [31,86, 120, 170,201, 21 4, 21 6]. Also left out 
is the vast area of numerical methods of solving differential equations, making 
use of their symmetry properties [80] or the treatment of asymptotic symmetries 
for difference equations [76] . A very active area of research, not covered in the 
present review, is the use of higher symmetries to identify integrable lattice 
equations [12, 13, 158-161, 223, 240] . 

Since we are reviewing a relatively new area of research, many open ques- 
tions remain. Thus, for purely difference equations on fixed lattices, the role of 
discrete symmetries has not been fully explored. Basically, what is needed, spe- 
cially for partial difference equations, is to apply results from crystallography, 
to characterize discrete or finite transformations taking a lattice into itself. For 
differential-difference equations with three independent variables, one or two of 
them discrete, a classification of equations with Kac-Moody-Virasoro symme- 
tries should help to identify new integrable lattice equations. Applications to 
genuine physical systems would be of great interest. The umbral approach of 
Section |21 has so far been applied in a rather formal manner. The question of 
the convergence of formal power series solutions must be addressed. The main 
question is whether one can develop a complete convergent quantum theory on 
a lattice. The generalized symmetries of Section 0] are an integral part of the 
theory of integrability on a lattice. There the greatest challenge lies in the field 
of applications, i.e in applying the techniques of integrability to the real world of 
discrete phenomena. Finally, the greatest challenge in the direction of symmetry 
adapted discretizations is to establish whether they provide improved numerical 
methods, in particular for partial differential equations, or higher order ordinary 
ones. 
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